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1.0  SUMMARY 


Experimental,  theoretical  and  numerical  investigations  of  quantum  computation  using  photon- 
based  qubits  (Quantum  Bits)  were  conducted  to  explore  the  Cluster  State  (or  one-way)  Quantum 
Computing  paradigm.  This  report  describes  research  on  a  unique  type  II  SPDC  (Spontaneous 
Parametric  Down  Conversion)  source  (“Schioedtei”)  design  that  can  generate  up  to  six  pairs  of 
entangled  photons  per  pass  through  the  type  II  crystal  assembly.  This  source  is  currently  being 
used  as  the  entangled  photon  source  to  create  photon-based  qubit  cluster  states.  Under  this 
project  we  developed  a  new  detector  design  architecture  that  turns  the  single  photon  detector  into 
a  number-resolving  detector  by  means  of  a  novel  three-dimensional  architecture  that  utilizes 
spatial  multiplexing.  We  have  studied  the  CNOT  gate,  as  an  archetypical  quantum  linear  optical 
gate,  and  found  several  interesting  features  in  the  both  the  ideal  and  the  realistic  case  of 
implementation  with  imperfect  (non-unit)  fidelity.  We  have  conducted  a  theoretical  investigation 
of  the  limitations  of  quantum  correlations  under  the  physically  imposed  constraint  of  no¬ 
signaling  (i.e.  no  faster  than  light  communication).  Finally,  we  discuss  our  ongoing  work  of 
quantum  algorithm  development  utilizing  the  measurement-based  quantum  computation 
approach,  and  compare  and  contrast  it  with  the  standard  quantum  circuit  model  approach  in  the 
case  of  an  unstructured  search. 

2.0  INTRODUCTION 

Under  this  AFRL/RI  in-house  project  we  continued  research  and  development  of  a  novel  multi¬ 
qubit  entangled  photon  source,  begun  under  the  AFRL/RI  in-house  project  “Quantum 
Information  Science  (QIS),”  (AFRL-RI-RS-TR-2012-073),  and  begun  investigations  into  the 
measurement-based  cluster  state  quantum  computation  paradigm  utilizing  photon-based  qubits. 
These  investigations  included:  (i)  the  development  and  characterization  of  a  new  multipli- 
entangled  photon  source  that  increased  the  usable  number  of  photon  pairs  by  a  factor  of  six  over 
conventional  entangled  photon  sources,  (ii)  design  of  multi-layer  superconducting  number¬ 
resolving  photon  detector,  (iii)  a  theoretical  and  experimental  investigation  into  the  requirements 
of  imperfect  (non-unit  fidelity)  two-qubit  linear  optical  photonic  gates,  and  (iv)  a  theoretical 
investigation  of  entanglement  and  nonlocality  addressing  the  issue  of  why  nature  does  not  take 
advantage  of  the  algebraically  allowed  maximum  correlations  amongst  collections  of  qubits.  In 
addition,  this  report  discusses  upgrades  to  our  in-house  AFRL/RI  Quantum  Computing 
Laboratory  under  the  current  project  and  our  thrust  to  transition  our  development  of  quantum 
gates/circuits  in  bulk  optics  to  an  on  chip  integrated  waveguide  implementation.  Finally,  this 
report  briefly  discusses  research  we  initiated,  and  is  currently  ongoing,  in  the  area  of  cluster  state 
quantum  algorithm  development. 

Cluster  State  Quantum  Computation  Background 

In  the  standard  Quantum  Circuit  Model  (QCM)  paradigm,  quantum  computations  are  executed 
by  successive  unitary  operations  acting  upon  an  initial  quantum  state  composed  of  many  qubits. 
These  unitary  operators  create  entanglement  amongst  the  qubits  through  quantum  interference. 
Entanglement  is  a  uniquely  non-classical  property  of  quantum  mechanical  systems  in  which  the 
correlations  between  sub-systems  can  be  stronger  than  that  allowed  by  classical  (conventional) 
computing  systems.  Recently  a  new  alternative  paradigm  for  quantum  computation  has  emerged 
called  One-way  Quantum  Computation  (OWQC)  [RuassendorfOl].  In  the  one-way  quantum 
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computer,  information  is  processed  by  sequences  of  single-qubit  measurements.  These 
measurements  are  performed  on  a  universal  resource  state — the  2D-cluster  state — which  does 
not  depend  on  the  algorithm  to  be  implemented.  The  new  approach  to  quantum  computation  goes 
by  the  collective  name  measurement-based  quantum  computation  (MBQC)  [Briegel09],  The 
appeal  of  MBQC  is  that  deterministic  quantum  computation  is  possible  based  on  (i)  the 
preparation  of  an  initial  entangled  cluster  state,  followed  by  (ii)  a  temporally  ordered  pattern  of 
single  qubit  measurements  and  feed-forward  operations  which  depend  on  the  outcome  of  the 
previously  measured  qubits  [RaussendorfOl],  Our  interests  in  OWQC  is  in  the  utilization  of 
photon-based  cluster  states  as  gates  and  circuits  for  quantum  computation  (see  [Vallone08],  and 
references  therein).  It  has  been  claimed  that  the  use  of  cluster  states  can  substantially  reduce  the 
resource  overhead  in  the  standard  QCM  to  photon-based  quantum  computation. 

In  the  OWQC  approach  a  quantum  computation  proceeds  as  follows:  (i)  A  classical  input  is 
provided  which  specifies  the  data  and  the  program,  (ii)  A  2D-cluster  state  of  sufficiently  large 
size  is  prepared.  The  cluster  state  serves  as  the  resource  for  the  computation,  (iii)  A  sequence  of 
adaptive  one-qubit  measurements  is  implemented  on  certain  qubits  in  the  cluster.  In  each  step  of 
the  computation  the  measurement  bases  depend  on  the  specific  program  under  execution  and  on 
the  outcomes  of  previous  measurements.  A  simple  classical  computer  is  used  to  compute  which 
measurement  directions  have  to  be  chosen  in  every  step,  (iv)  After  the  measurements  the  state  of 
the  system  has  the  product  form  where  a  indexes  the  collection  of  measurement 

outcomes  of  the  different  branches  of  the  computation.  The  states  |  t//";  ^  in  all  branches  are  equal 

to  the  desired  output  state  up  to  a  local  (Pauli)  operation.  The  measured  qubits  are  in  a  product 
state  1 4^  which  also  depends  on  the  measurement  outcomes.  The  OWQC  is  computationally 

universal,  i.e.  even  though  the  results  of  the  measurements  in  every  step  of  the  computation  are 
random,  any  quantum  computation  can  deterministically  be  realized.  Notice  that  the  temporal 
ordering  of  the  measurements  plays  an  important  role  and  has  been  formalized  as  a  feed-forward 
procedure  [RaussendorfOl], 

In  realistic  physical  systems  decoherence  tends  to  make  quantum  systems  behave  more 
classically.  One  could  therefore  expect  that  decoherence  would  threaten  any  computational 
advantage  possessed  by  a  quantum  computer.  However,  the  effects  of  decoherence  can  be 
counteracted  by  quantum  error  correction  [Shor96].  In  fact,  arbitrarily  large  quantum 
computations  can  be  performed  with  arbitrary  accuracy  provided  the  error  level  of  the 
elementary  components  of  the  quantum  computer  is  below  a  certain  threshold.  This  important 
result  is  called  the  threshold  theorem  of  quantum  computation  [Aliferis06]. 

Fault-tolerant  schemes  for  OWQC  using  photons  have  recently  been  developed  [Dawson06, 
Vamava06].  The  dominant  sources  of  error  in  this  setting  are  photon  loss  and  gate  inaccuracies. 
The  constraint  of  short-range  interaction  and  arrangement  of  qubits  in  a  2D  lattice — a 
characteristic  feature  of  the  initial  one-way  quantum  computer — is  not  relevant  for  photons.  In 
[Dawson06]  both  photon  loss  and  gate  inaccuracies  were  taken  into  account  yielding  a  trade-off 
curve  between  the  two  respective  thresholds.  Fault- tolerant  optical  computation  is  possible  for  a 
gate  error  rate  of  1(T4  and  photon  loss  rate  of  3xl0"3.  In  [Vamava06]  the  stability  against  the  main 
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error  source  of  photon  loss  was  discussed.  With  non-unit  efficiencies  r|S  and  r|D  of  photon 
creation  and  detection  being  the  only  imperfections,  the  very  high  threshold  of  r|sr|D  >  2/3  was 
established.  Further,  encoding  a  collection  of  physical  qubits  within  the  2D  cluster  state  offers  a 
means  of  topological  error  protection  for  the  logical  qubit.  Topologically  protected  quantum  gate 
operations  are  performed  by  measuring  some  regions  of  qubits  in  the  Z-basis,  which  effectively 
removes  the  qubits  from  the  state.  The  remaining  cluster,  whose  qubits  are  measured  in  the  X  - 
and  X  ±Y  -  bases,  thereby  attains  a  non-trivial  topology  in  which  fault-tolerant  quantum  gates 
can  be  encoded.  A  topological  method  of  fault-tolerance  for  OWQC  can  then  be  achieved 
[Raussendorf07]. 

Experimental  Research: 

Photons  are  particularly  desirable  for  quantum  information  processing  tasks  since  they  are 
relatively  free  from  environmental  decoherence.  Hence,  they  are  also  essential  for  any  long 
distance  conveyance  of  quantum  information,  and  do  not  require  cryogenic  cooling.  Entangled 
photon  sources  with  the  highest  mode  quality  are  based  on  spontaneous  parametric  down 
conversion  (SPDC).  This  is  a  process  where  laser  pump  photons  are  converted  into  ‘signal’  and 
‘idler’  entangled  pairs  in  nonlinear  (NL)  crystals.  SPDC  in  nonlinear  crystals  has  provided  the 
optical  sources  for  groundbreaking  foundational  and  applications  work  in  quantum  optics  (QO) 
for  the  last  two  decades  [O’Brien07], 

SPDC  is  an  inherently  inefficient  process,  and  work  based  on  it  is  generally  limited  by  the  net 
signal  level  or  the  number  of  photons  that  can  be  entangled  in  given  applications.  Photon  yield  is 
related  to  laser  power,  which  cannot  be  increased  beyond  the  level  where  higher  order  NL 
contributions  (multi-photon  events)  yield  errors  in  quantum  processing  applications.  This  point 
has  now  been  reached  in  applications  that  require  independent  sources  of  entangled  qubits.  The 
work  begun  under  the  in-house  Quantum  Information  Science  project  focused  on  (i)  developing 
a  6-qubit  capable  photon-based  quantum  information  testbed  and  (ii)  initial  development  of  new 
sources  of  entangled  photons  that  greatly  increase  process  efficiency,  without  increasing  laser 
power,  in  a  regime  where  high  detection  quantum  efficiency  is  available  -  a  highly  desirable  goal 
not  previously  accomplished  in  the  scientific  community  to  date.  This  latter  direction  of  research 
was  continued  and  expanded  upon  in  the  current  in-house  project  Cluster  State  Quantum 
Computing. 

Number  resolving  photon  counting  at  the  single  photon  level,  i.e.  distinguishing  1,  2  or  3  photons 
is  an  important  experimental  ability.  While  experiments  can  be  performed  without  number 
resolving  detectors  such  as  the  APD  (Avalanche  Photo  Diode)  we  are  currently  using,  a 
significant  number  of  interesting  experiments  require  the  number  resolving  ability.  Therefore  we 
considered  a  known  single  photon  detector  (click  or  no-click)  and  have  developed  a  new  detector 
design  architecture  that  turns  the  single  photon  detector  into  a  number-resolving  detector.  This  is 
done  through  spatial  multiplexing.  The  architecture  we  have  developed  allows  for  higher  density 
of  detection  elements  and  larger  number  of  detector  elements  then  the  current  state  of  the  art. 
This  leads  to  faster  repetition  times,  and  most  importantly,  superior  number  resolution. 


Approved  for  Public  Release;  Distribution  Unlimited. 

3 


Theory/Numerical  Research: 

The  creation  of  single  photons  is  most  often  performed  by  non-linear  processes  such  as  SPDC. 
However  to  create  cluster  states  from  these  resource  photons  requires  linear  optics,  such  as  the 
CNOT  or  CZ  gates.  Numerous  implementation  methods  have  been  suggested  for  these  gates  in 
theory,  but  in  practice  any  implementation  will  be  imperfect.  Therefore  it  is  important  to  be  able 
to  characterize  and  optimize  imperfect,  i.e.  fidelity  less  than  1,  linear  optical  gates  and  state 
transformations.  We  have  studied  the  CNOT  gate,  as  an  archetypical  linear  optical  gate  and 
found  several  interesting  features  in  the  both  the  ideal  case  and  the  realistic  case  of  imperfect 
fidelity.  In  particular  we  find  that  the  success  rate  can  be  increased  as  the  fidelity  decreases. 
This  trade  off  in  fidelity  for  success  opens  several  interesting  possibilities  in  the  field  of  linear 
optics.  We  also  found  that  the  CNOT  gate  has  a  high  degree  of  symmetry  that  simplifies  its 
physical  implementation  in  both  cases  and  proposed  a  practical  experiment  to  test  our  theories. 
Such  studies  can  and  have  been  carried  over  to  other  gates  and  state  transformation  with  relative 
ease. 

The  ultimate  goal  of  MBQC  is  to  execute  quantum  algorithms  with  a  speedup  over  their 
corresponding  classical  algorithm  counterparts.  Under  this  project  we  began  an  investigation  of 
Grover’s  search  algorithm  (GSA)  on  an  unsorted  list  of  elements  in  the  MBQC  paradigm,  and  its 
comparison/contrast  with  the  usual  QCM  approach.  GSA  serves  as  an  important  prototypical 
benchmark  for  many  numerical  simulations  of  quantum  algorithms  [Grover97,  Walther05].  Our 
preliminary  results  (which  are  currently  being  written  up  for  journal  submission)  indicate  that  the 
MBQC  implementation  of  Grover’s  algorithm  is  faster  than  even  Grover’s  quantum  algorithm 
(with  its  quadratic  speedup  over  a  brute  force  search). 

In  brief,  Grover’s  oracle-based  unstructured  search  algorithm  is  often  stated  as  “given  a  phone 
number  in  a  directory,  find  the  associated  name.”  More  formally,  the  problem  can  be  stated  as 
“given  as  input  a  unitary  black  box  Uj  for  computing  an  unknown  function  /:  [0,1 }"  — > {0,1  [find 
x=xo  an  element  of  {0,1}"  such  that  /(xo)  =1,  (and  zero  otherwise).”  The  crucial  role  of  the 
externally  supplied  oracle  Uf  (whose  inner  workings  are  unknown  to  the  user)  is  to  change  the 
sign  of  the  solution  |*0),  while  leaving  all  other  states  unaltered.  Thus,  Ut  depends  on  the  desired 

solution  xo.  Under  the  previous  in-house  QIS  project,  we  developed/simulated  an  amplitude 
amplification  algorithm  in  which  the  user  encodes  the  directory  (e.g.  names  and  telephone 
numbers)  into  an  entangled  database  state,  which  at  a  later  time  can  be  queried  on  one  supplied 
component  entry  (e.g.  a  given  phone  number  to)  to  find  the  other  associated  unknown  component 
(e.g.  name  xo).  For  N=2"  names  |x)  with  N  associated  phone  numbers]/),  performing  amplitude 

amplification  on  a  subspace  of  size  N  of  the  total  space  of  size  N2  produces  the  desired  state 

|x0)|/0)in  VjV  steps. 

3.0  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 

3.1.  Multipli-entangled  photons  from  a  spontaneous  parametric  down-conversion  source 

Photon  based  quantum  computation,  with  single  or  entangled  photons,  is  a  heavily  researched 
area.  This  is  in  part  due  to  many  desirable  properties  of  photons  such  as  (i)  room  temperature 
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operation,  (ii)  immunity  from  the  environment  and  (iii)  superior  mode  quality.  Spontaneous 
parametric  down  conversion  has  proven  to  be  the  most  reliable  method  of  generating  entangled 
photon  pairs.  Type  I  sources  spontaneously  convert  one  linearly  polarized  parent  photon  into  two 
daughters,  each  having  a  polarization  orthogonal  to  the  parent.  The  spontaneous  nature  of 
parametric  down  conversion  produces  a  ring  pattern  where  each  diametric  photon  pair  shares  the 
same  parent  photon.  The  type  I  process  produces  two  photons  of  the  same  polarization  which  are 
path  entangled.  That  is,  detecting  a  photon  (signal)  in  path  A  implies  that  its  sister  (idler)  can  be 
found  in  the  diametrically  opposite  spot  [DragomanOl].  This  implies  that  in  the  polarization 
basis  a  mixed  state  \H)\\H)2  or  |T)i |V)2  will  be  produced.  Many  experiments  that  require  photon 
pairs,  but  not  entangled  pairs,  use  the  output  of  a  type  I  crystal  as  input  to  a  more  sophisticated 
experiment.  Type  I  crystals  are  inherently  birefringent,  but  the  walk-off  associated  with  these 
crystals  is  mitigated  by  the  fact  that  the  down  converted  photons  are  the  same  polarization;  the 
delays  that  the  signal  and  idler  photons  experience  are  the  same.  Type  I  sources  have  been  used 
for  many  years  in  harmonic  generation  (SHG,  THG)  systems  as  frequency  converters. 

Kwiat  first  described  a  feasible  source  for  SPDC-generated  entangled  pairs  using  type  I  P-BBO 
[Kwiat99]  (beta-Barium  borate,  BaB204).  This  consisted  of  a  stacked  pair  of  type  I  crystals 
rotated  90°  relative  to  each  other  (Figure  1).  This  allows  for  the  generation  of  two  orthogonally 
polarized  cones  that  overlap  in  space.  Each  crystal  can  only  be  excited  by  a  certain  linear 
polarization.  The  stack  must  be  pumped  by  a  beam  made  up  of  components  that  excite  each 
crystal  equally.  That  is,  if  the  stack  consists  of  an  optic  axis  that  is  vertical  in  the  first  crystal,  and 
horizontal  in  the  second,  the  pump  beam  polarization  must  be  oriented  at  45°.  The  resulting 
superposition  state  is  |//)i |//)2  ±ew  |l/)i  |V)2.  It  is  important  to  note  that  there  is  a  temporal  delay 
associated  with  the  down  converted  states  due  to  the  crystal’s  birefringence.  Compensation 
depends  on  the  wave  packet  of  the  pump  photon. 


Figure  1.  Kwiat’ s  type  I  pair.  Each  crystal  will  spontaneously  down  convert  a  linearly  polarized  photon 
pair  orthogonal  to  its  pump  photon.  When  specifically  oriented,  these  down  converted  rings  can  overlap 
and  create  a  polarization  entangled  pair. 

Type  II  sources  [Kwiat95]  down  convert  a  linearly  polarized  parent  photon  into  two  orthogonally 
polarized  daughters  making  them  particularly  interesting  because  the  crystal  is  birefringent.  One 
daughter  photon  will  walk  off  faster  than  the  other  and  lead  to  a  noticeable  spatial  separation,  and 
thus  there  are  two  intersecting  cones  (Figure  2).  The  walk  off  of  type  II  crystals  limits  the  length 
of  the  crystal  because  the  extraordinary  index  of  refraction  will  quickly  bend  light  out  of  the 
crystal.  Indistinguishable  photons  are  produced  in  the  intersections  of  the  two  cones.  These  two 
points  form  a  superposition  state  of  polarization  (|//)i|F)2  ±  ew  |k)i|//)2)  and  have  been 
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exhaustively  studied  and  used  as  inputs  to  more  complex  photonic  systems.  Figure  3  shows  the 
evolution  of  the  rings  of  entangled  photons  produced  from  SPDC  crystals  under  type  I  and  II 
phase  matching  conditions. 


Figure  2.  Standard  type  II  down  conversion.  A  linear  pump  beam  spontaneously  down  converts  to  two 
photons,  one  of  which  has  the  same  polarization  as  the  pump.  The  other  is  orthogonal.  The  familiar 
double  ring  pattern  is  a  product  of  the  crystal’s  birefringence. 


Figure  3:  Imaging  of  the  down-converted  light  for  three  different  configurations.  First  row:  type-I 
SPDC  as  a  function  of  the  tilt  of  one  crystal.  Second  row:  type-I  SPDC  rings  of  different  diameters 
as  a  function  of  the  polarization  of  the  pump  beam  (horizontal  on  the  left  and  vertical  on  the  right). 

Third  row:  type-II  SPDC  rings  as  a  function  of  the  tilt  of  the  crystal.  All  cases  involve  the  CW 
pump  laser  beam. 

In  a  similar  fashion  to  Kwiat’s  type  I  stack,  Bitton  et  al.  [BittonOl]  described  a  type  II  stack 
comprised  of  two  crystals  rotated  180°  relative  to  each  other  (Figure  4).  This  allows  the  linear 
pump  scheme  to  remain  unchanged  and  yields  one  set  of  rings  from  either  crystal.  The  set  of 
rings  entirely  overlap  each  other  and  thus  can  yield  an  entangled  photon  pair  of  the  same  state  as 
standard  type  II.  Addressing  the  compensation  is  a  necessary  requirement  with  any  birefringent 
crystals.  A  standard  type  II  stacked  configuration  allows  for  greater  pair  production  and  more 
useable  detection  area.  In  this  source  as  well  as  a  type  I  stack,  the  fundamental  size  of  the 
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collection  apertures  become  the  limiting  factor  in  the  number  of  entangled  pairs  that  can  be 
collected. 


Figure  4.  Type  II  crystal  assembly  as  described  by  [BittionOl],  The  second  crystal  is  rotated  180°  relative 
to  the  first,  resulting  in  two  overlapping  sets  of  cones  with  orthogonal  polarization. 

U’ren  et  al.  [U’ren06]  described  a  type  II  crystal  assembly  (Figure  5a)  that  is  designed  for  group 
velocity  matching  (GVM)  of  the  pump  and  signal/idler  wave  packets,  thereby  removing  any 
spectral  distinguishability  of  the  down  converted  photons.  The  assembly  consists  of  a  successive 
stack  of  nonlinear  crystals  (P-BBO,  BiBO  (Bismuth  Borate,  BiB306))  separated  by  a  thin  layer  of 
compensating  crystal  (calcite  (CaC03),  a-BBO).  This  slowly  compensates  different  components 
(pump  and  down  converted  wave  packets)  such  that  by  the  end  of  the  stack  there  is  no  spectral 
walk  off.  The  need  to  spectrally  filter  post  down  conversion  is  mitigated  by  the  symmetry  of  their 
joint  spectral  function  (Figure  5b).  Removing  this  requirement  typically  increases  the  useable 
count  rate  and  overall  efficiency. 


Figure  5.  (a)  Type  II  custom  assembly  showing  alternating  layers  of  P-BBO  (red)  and  calcite  (blue),  (b) 
Joint  spectral  function  of  down  converted  wavepacket.  This  implies  a  maximally  separable  state. 


Type  I  and  II  crystals  are  still  governed  by  their  spontaneous  nature,  and  this  becomes 
problematic  when  large  numbers  of  entangled  photons  are  required.  In  a  typical  configuration  for 
the  generation  of  greater  than  four  photons  a  cascaded  apparatus  is  used.  For  this  setup  either 
multiple  crystals  are  used  in  succession,  or  multiple  passes  through  a  single  crystal  (Figure  6). 
This  implies  an  overall  increase  in  footprint  size.  Hyper-entanglement  has  been  considered  to 
mitigate  the  spontaneous  nature  of  down  conversion  by  adding  entanglement  various  degrees  of 
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freedom,  not  just  polarization  [Ceccarelli09].  While  this  is  effective  it  requires  larger  physical 
hardware  requirements  and  more  complicated  analysis  processes. 


Figure  6.  Cascaded  and  multi-pass  crystal  configurations  for  the  generation  of  cluster  states  [Lu07], 

In  recent  years  there  has  been  a  paradigm  shift  in  quantum  computation  with  the  need  to  migrate 
toward  schemes  that  require  only  single  qubit  measurements.  One-way  quantum  computation 
(cluster  state)  has  facilitated  this  shift.  Cluster  state  computation  allows  a  predetermined 
sequence  of  single  qubit  measurements  to  determine  the  algorithm  being  evaluated  [Walther05], 
This  protocol  requires  a  highly  entangled  cluster  state  [RaussendorfOl]  generated  from  a 
resource  of  qubits.  Such  a  cluster  state  can  be  constructed  by  preparing  each  of  the  qubits  into  a 
state,  |  +)  =  T=(|0)  +|1)),  and  applying  controlled-phase  gates  to  link  the  required  qubits. 
Computation  proceeds  with  a  sequence  of  single  qubit  measurements  whose  results  will 
classically  feed  forward  to  control  the  basis  required  for  future  measurements  [Nielson05]. 
Cluster  state  computation  allows  for  a  practical  resource  reduction  in  qubits  and  hardware 
compared  to  other  quantum  computing  methods.  That  being  said,  the  fundamental  requirement 
for  larger  numbers  of  qubits  still  exists.  The  source  we  have  developed  produces  larger  qubits 
numbers  than  that  of  a  typical  type  II  SPDC  source. 

SPDC  custom  crystal  assembly 

Our  custom  two-crystal  assembly  (designated  as  “Schioedtei”  henceforth)  design  consists  of  a 
pair  of  type  II  non-collinear  phase-matched  SPDC  crystals  cut  for  degenerate  down-conversion 
whose  optic  axes  are  rotated  orthogonal  with  respect  to  one  another.  The  pair  of  crystals  is 
optically  contacted  with  one  another  and  a  dual  band  (405/810  nm)  anti-reflection  coating 
applied  to  the  two  exterior  faces  of  the  assembly.  Any  type  II  material  can  be  used  to  create  an 
equivalent  device.  Our  particular  version  that  will  be  discussed  here  was  constructed  from  two 
8x8x2  mm  type  II  beta-Barium  borate  (P-BBO,  BaB204)  crystals  phase  matched  (at  angles  of 
theta  =  41.9°,  phi  =  30°)  for  810  nm  spontaneous  parametric  down-conversion. 

Exciting  Schioedtei  with  an  incident  45°  polarized  pump  beam  produces  one  pair  of  rings  from 
each  of  the  type  II  crystals.  Each  pair  of  rings  is  orthogonal  to  the  other  resulting  in  12 
intersection  points  (or  simply  “points”)  where  indistinguishable  photons  are  produced.  Referring 
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to  Figure  7,  the  indicated  points  marked  5,  6  (Bell  pair  #1  from  crystal  #1)  and  7,  8  (Bell  pair  #2 
from  crystal  #2)  are  the  typical  Bell  states,  \ip)5, 6(7, 8) =  ( | H V)5, 6(7, 8)  ±  | The 


Figure  7.  Type-II  SPDC  Schioedtei  source.  See  text  for  discussion  of  the  intersection 
points  of  the  overlapping  rings. 

points  indicated  by  1,  2,  3,  4  are  the  product  of  two  bell  states,  \ip)i,2,3A  =  ^  (|HF)i,4  ± 

|]///)L4)  (\HV)23  ±  e~l<p\VH)2.i),  produced  from  photons  from  both  crystal  1  and  2 
concurrently.  Points  9,  11  and  10,  12  are  |PV)9,n  and  |//H)io,i2  states  produced  from  photons 
from  crystal  1  and  2  concurrently.  Further  analysis  and  experimental  results  of  Schioedtei  are 
covered  in  section  4. 1 . 

3.2  A  multi-layer  three  dimensional  superconducting  nanowire  photon  detector 

Construction  of  photon-counting  devices  with  high  counting  efficiency,  high  number  resolution 
and  short  reset  times,  is  highly  desirable  for  a  wide  array  of  applications,  such  as  quantum  key 
distribution  [Xu08],  quantum  communication  [O’Brein09],  quantum  computing  [KnillOl], 
[UskovlO],  [Knill02]  among  others  [Hadfield09],  [Smithll],  Here  we  describe  and  perform 
some  simple  analyses  of  a  proposed  detector  design  that  uses  multiple  short  sections  of 
superconducting  nanowires  to  construct  a  new  superconducting  nanowire  single  photon  detector 
(SNSPD).  We  refer  to  these  short  sections  of  nanowire  as  pixels  and  arrange  them  in  a  two 
dimensional  grid  in  analogy  with  a  standard  CCD  camera.  We  will  discuss  the  potential 
advantages  of  such  a  system  and  the  difficulties  of  the  design. 

When  an  incident  photon  strikes  a  Niobium  nitride  (NbN)  nanowire,  or  other  superconducting 
material  such  as  NbTiN  or  a-WxSii_x  developed  recently  at  NIST,  it  creates  a  resistive  hot  spot 
[Namll],  This  hot  spot  causes  the  current  in  the  superconductor  to  deflect  around  the  spot,  thus 
increasing  the  current  density  in  the  wire.  This  increased  current  density  leads  to  an  increase  in 
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the  temperature  of  a  small  section  of  the  wire.  If  the  nanowire  is  held  just  below  the  critical 
current  for  superconduction,  then  the  increase  in  heat  will  break  the  superconducting  condition 
and  the  resistance  of  the  wire  will  spike  upward  for  a  short  time.  This  resistance  spike  creates  a 
measurable  current  in  the  external  resistance  load  and  a  photon  is  counted. 

Present  superconducting  nanowire  systems,  such  as  NbN,  have  reasonably  good  counting 
efficiency  [DaulerlO],  [Marsilill],  by  which  we  mean  the  probability  of  an  incident  photon 
being  detected  is  over  25%.  However,  a  significant  problem  exists  with  the  number  resolution, 
relaxation  time,  and  fill  factors  [Gurevich87],  [DaulerlO],  [Marsilill],  A  detector  consisting  of  a 
single  wire  can  be  made  to  cover  a  significant  detection  area  by  creating  a  meander.  Usually  this 
means  folding  the  nanowire  back  and  forth  across  the  desired  area  of  approximately  10  pm  x  10 
pm  [DaulerlO],  [Marsilill],  This  however  is  not  a  number  resolving  detector.  All  that  the 
detector  can  feel  is  the  loss  of  the  superconducting  condition  somewhere  in  the  nanowire.  Should 
two  photons  strike  the  wire  simultaneously  in  two  different  locations  the  current  drop  is  very 
similar.  One  suggested  solution  to  this  lack  of  number  resolving  capability  is  to  increase  the 
number  of  wires  in  the  meander.  This  has  been  done  experimentally  by  Dauler  et  al  [DaulerlO], 
While  this  approach  improves  on  the  single  wire  meander  it  still  consists  of  long  wires  each  of 
which  occupies  a  significant  portion  of  the  active  detection  area  (i.e.  each  nanowire  in  a  4 
nanowire  meander  takes  approximately  25%  of  the  active  area).  In  order  to  have  a  high 
probability  of  correctly  detecting  n  number  of  photons  one  would  need  significantly  more  than  n 
wires. 

We  proposed  a  detector  design  using  short  sections  of  wire,  which  we  will  refer  to  as  pixels,  that 
are  arranged  in  a  2D  grid  to  create  the  detection  area.  Such  a  design  would  use  a  large  number  of 
pixels  thus  giving  high  number  resolution  and  the  small  size  of  the  pixels  gives  short  relaxation 
times.  We  call  this  configuration  a  multi-layer  superconducting  number-resolving  photon 
detector.  A  significant  problem  with  creating  a  two  dimensional  array  of  nanowire  pixels  is  the 
question  of  how  to  attach  the  leads  to  each  pixel,  as  the  leads  are  of  a  similar  size  as  the  pixels 
themselves.  One  could  simply  move  the  pixels  farther  and  farther  apart  to  fit  in  all  the  necessary 
connections  but  this  is  impractical  as  the  space  between  pixels  does  not  detect  photons  and  the 
device’s  overall  efficiency  would  decrease  below  useful  levels.  Ideally  the  pixels  will  be  packed 
as  closely  as  possible,  while  still  avoiding  cross  talk.  This  will  maximize  the  so  called  fill  factor, 
the  ratio  of  the  photon  sensitive  area  to  the  non-sensitive  area  within  the  “active”  area  of  the 
detector.  We  therefore  propose  moving  away  from  the  two  dimensional  approaches  used  to  date 
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Figure  8.  A  single  superconducting  nanowire  “pixel”  bridge. 


~5nm 


and  instead  suggest  a  multilayer,  three  dimensional  architecture.  In  this  scheme  the  non¬ 
superconducting  leads  are  allowed  to  pass  under  the  active  detector  pixels.  To  create  this  effect 
we  shape  the  pixels  like  small  bridges  as  seen  in  Figure  8.  This  shape  was  chosen  because  of  its 
relatively  simple  design  and  in  order  to  maximize  the  fill  factor. 

This  creates  a  three  layer  design,  with  the  bottom  layer  containing  the  leads,  an  insulating  middle 
layer,  and  the  active  detection  layer  on  top.  We  now  show  a  bird’s  eye  view  (plan  view)  of  the 
final  device  with  all  three  layers  aligned  on  top  of  each  other,  Figure  9.  The  black  arrows  show 
the  movement  of  the  current  throughout  the  device. 


Figure  9.  Plan  view  of  final  device  with  all  three  layers  aligned  on  top  of  each  other. 
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The  current  enters  the  detector  in  the  bottom  layer  (red)  through  the  shared  input  lead  (center  of 
figure).  Current  then  moves  up  through  the  middle  layer  connections,  called  “posts"  (green),  to 
the  top/detection  layer.  Once  in  the  top  layer  it  moves  along  the  surface  of  the  bridge  (purple). 
This  is  the  area  in  which  an  incident  photon  will  form  a  resistance  blockage.  The  current  then 
moves  back  down  to  the  bottom  layer  and  is  channeled  out  of  the  device  by  the  output  leads 
(red).  Note  that  the  leads  (red)  pass  under  the  pixel  bridges,  between  the  posts  and  that  the 
input/output  leads  are  all  on  the  same  layer.  The  external  detection  electronics  would  then  be 
similar  industry  standards.  As  a  final  aside  each  pixel  can  be  wired  as  a  completely  independent 
circuit,  but  the  number  of  leads  will  increase  and  counting  simultaneous  events  between  elements 
can  become  difficult. 

3.3  Theory/experimental  requirements  of  imperfect  two-qubit  linear  optical  photonic  gates 

Knill,  Laflamme,  and  Milbum  (KLM)  significantly  advanced  the  prospect  of  single-photon 
quantum  computing  in  their  seminal  paper  [KnillOl],  in  which  they  overcame  the  need  for 
nonlinear  interactions  by  using  the  inherent  nonlinearity  of  photon  measurements.  In  this 
scheme,  the  computational  system  is  combined  with  ancillary  modes,  and  the  gate  operation  is 
performed  on  the  enlarged  state  space.  The  ancilla  modes  are  measured  with  photon-number- 
resolving  detectors,  such  as  those  described  above  leaving  the  computational  modes  undisturbed 
and  in  the  desired  output  state  provided  the  measurement  is  successful.  In  our  previous  work 
[UskovlO,  Uskov09,  Smithll],  we  have  shown  that  a  combination  of  analytical  and  numerical 
techniques  may  be  used  to  design  optimal  linear  optical  transformations  implementing  two-  and 
three-qubit  entangling  gates.  Here  we  show  results  for  non-ideal  gates  and  suggest  an  experiment 
to  test  them. 

The  probabilistic  nature  of  quantum  measurement  implies  a  trade-off  between  the  success  rate  of 
the  operation  (the  probability  of  obtaining  the  desired  measurement  outcome  for  the  ancillary 
modes)  and  the  fidelity  (the  overlap  between  the  actual  and  desired  states  of  the  computational 
system  when  the  ancilla  measurement  is  successful).  Previously,  solutions  were  obtained  that 
have  the  maximum  possible  ancilla  measurement  success  probability  given  the  constraint  of 
perfect  fidelity  for  a  specified  transformation  [UskovlO,  Uskov09],  In  practical  implementations, 
however,  the  goal  of  perfect  fidelity  may  not  always  be  desirable  or  even  attainable.  We  have 
therefore  generalized  our  previous  techniques  to  the  case  of  imperfect  fidelity,  and  investigated 
the  above-mentioned  trade-off  between  the  fidelity  and  success  of  the  linear  optical 
transformations.  It  was  found  that  for  sufficiently  small  deviations  from  perfect  fidelity,  a  single 
optimization  parameter  determines  the  relationship  between  fidelity  and  optimal  success  rate 
[Smithll], 

The  input  state  to  the  experiment  |Tcomp>x|Tancilla>  is  a  product  of  the  computational  state 
containing  Me  photons  in  Nc  modes,  and  an  ancilla  state  containing  Ma  photons  in  Na  modes. 
The  Nc  computational  modes  are  those  on  which  the  actual  gate  is  intended  to  act.  Assuming 
dual-rail  encoding,  each  qubit  is  represented  by  one  and  only  one  photon  in  two  computational 
modes,  so  we  have  Me  =  Nc=2.  The  ancilla  state  may  in  general  be  separable,  entangled,  or  a 
maximally  entangled  or  ebit  state  carrying  spatially  distributed  entanglement  [Wilde09],  though 
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here  we  propose  using  only  a  product  state  of  single-photon  and  zero-photon  ancillas,  which  are 
relatively  simple  to  produce  in  an  experimental  setting. 

The  linear  optical  device  transforms  the  creation  operator  associated  with  each  input  mode 
i  to  a  sum  of  creation  operators  Y.j^i,jaj°Ut^ ■  Here  U,  which  contains  all  physical  properties  of 
the  device,  is  an  N  x  N  matrix,  where  N  =  Nc  +  Na  is  the  total  number  of  modes.  The  total  input 
state  may  be  written  as  a  superposition  of  Fock  states  \W  >=  \ nltn12, ... ,  nN  >,  where  n,  is  the 
occupation  number  of  the  i-th  input  mode,  and  Yjni  =  Mc  +  Ma  =  M  is  the  total  number  of 
photons.  The  input  state  is  transformed  as 

Fom  >  =  n\Vou,  >=  utJ  (1) 


We  note  that  Q  is  a  multivariate  polynomial  of  degree  M  in  the  elements  Uij.  Once  the 
transformation  is  complete,  a  measurement  is  applied  to  the  Na  ancillary  modes.  In  the  case  of  a 
number-resolving  photon-counting  measurement, (^measured.]  =  (^nc+ i> ^nc+ 2-  ■■■  -  KN |,  where  K, 
is  the  number  of  photons  measured  in  the  i-th  mode  of  the  ancilla.  The  resulting  transformation 
of  the  computational  state  is  a  contraction  quantum  map  =  ^|^0mP)/||^|^°mP)|| 

[Kraus83],  where  A  =  A(U)  is  defined  by, 

A^comp)  =  (Knc+1>Knc+2 . KN\n\Win).  (2) 


The  linear  operator  A,  which  maps  computational  input  states  to  computational  output  states, 
contains  all  the  information  of  relevance  to  the  transformation.  We  define  the  fidelity  as  the 
probability  that  the  desired  target  gate  ATar  has  been  faithfully  implemented  on  the  computational 
modes  given  a  successful  measurement  of  the  ancilla  modes: 


F(A)  = 


\Tr(AfATar)\2 
2McTr(A+d)  ’ 


(3) 


since  Tr  (ATari  ATar)=2Mc  for  a  properly  normalized  target  gate.  As  we  are  interested  in  deviations 
from  perfect  fidelity,  we  define  S  =1-F  as  our  main  parameter  [Smithll], 


We  define  the  success  rate  of  the  ancilla  measurement  to  be  given  by  an  average  over  all 
computational  input  states, 


,  ,  Tr(AUTar) 

5(4)  -  2Mc\\U\\2M  ’ 

for  general  complex  U.  Note  that  U  need  not  be  unitary,  as  any  matrix  can  be  made  unitary  via 
the  unitary  dilation  technique  by  adding  vacuum  modes  [Knill02,  uskov09].  We  also  note  that 
the  Hilbert-Schmidt  norm  {A\A)  =  Tr  (A^A)/2Mc  ,  used  in  our  definition  of  S,  is  bounded  above 
by  the  square  of  the  operator  norm,  ||i4||2  =  (||4||max)2  =  Max((lF^omp|dTartJ4|lf/^omp))  and 
below  by  (||A||mm)2  =  Min{J^^°mv\ATarA( A\W^°mv^)  where  the  maximum  and  minimum  are 
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taken  over  the  set  of  properly  normalized  input  states.  In  the  limit  F  — >1,  ||i4||?nax/||i4||?nm — »1, 
and  all  definitions  of  the  success  rate  coincide. 

3.4  Nonlocality,  entanglement  witnesses  and  supra-correlations 

While  entanglement  is  believed  to  underlie  the  power  of  quantum  computation  and 
communication,  it  is  not  generally  well  understood  for  multipartite  systems.  Recently,  it  has  been 
appreciated  that  there  exists  proper  no-signaling  probability  distributions  derivable  from 
operators  that  do  not  represent  valid  quantum  states.  Such  systems  exhibit  supra-correlations 
that  are  stronger  than  allowed  by  quantum  mechanics,  but  less  than  the  algebraically  allowed 
maximum  in  Bell-inequalities  (in  the  bipartite  case).  Some  of  these  probability  distributions  are 
derivable  from  an  entanglement  witness  W,  which  is  a  non-positive  Hermitian  operator 
constructed  such  that  its  expectation  value  with  a  separable  quantum  state  (positive  density 
matrix)  psep  is  non-negative  (so  that  Tr[  W  p]<  0  indicates  entanglement  in  quantum  state  p).  In 
the  bipartite  case,  it  is  known  that  by  a  modification  of  the  local  no-signaling  measurements  by 
spacelike  separated  parties  A  and  B,  the  supra-correlations  exhibited  by  any  W  can  be  modeled  as 
derivable  from  a  physically  realizable  quantum  state  p.  However,  this  result  does  not  generalize 
to  the  «-partite  case  for  ri> 2.  Supra-correlations  can  also  be  exhibited  in  2-  and  3-qubit  systems 
by  explicitly  constructing  “states”  O  (not  necessarily  positive  quantum  states)  that  exhibit  PR 
correlations  for  a  fixed,  but  arbitrary  number,  of  measurements  available  to  each  party.  In  this 
area  of  research  we  examined  the  structure  of  “states”  that  exhibit  supra-correlations.  In  addition, 
we  examined  the  affect  upon  the  distribution  of  the  correlations  amongst  the  parties  involved 
when  constraints  of  positivity  and  purity  are  imposed.  We  investigated  circumstances  in  which 
such  “states”  do  and  do  not  represent  valid  quantum  states. 

Physics  imposes  limits  on  the  correlations  that  can  be  observed  by  distant  (i.e.  spacelike 
separated)  parties.  In  particular,  special  relativity  (SR)  implies  the  principle  of  no-signaling 
(NS),  that  is,  correlations  cannot  lead  to  any  sort  of  instantaneous  communication  between 
spacelike  separated  observers.  Quantum  correlations  may  be  stronger  than  classical,  and  their 
violation  of  Bell  inequalities  (BI)  [Bell64]  suggests  that  quantum  mechanics  (QM)  cannot  be 
regarded  as  a  local  realism  theory.  Tsirelson  [Tsirelson80]  showed  that  there  is  an  upper  bound 
to  the  violation  of  BI,  which  implies  that  the  amount  of  non-locality  allowed  by  QM  is  limited. 
Popescu  and  Rohrlich  (PR)  showed  [Popescu94]  that  there  exists  a  broad  class  of  no-signaling 
theories  which  allow  stronger-than-quantum  or  supra-quantum  correlations.  PR  developed  a 
valid  joint  probability  distribution  whose  violations  of  the  BI  lie  above  those  of  physical 
quantum  correlations  and  below  the  allowed  algebraic  maximum  of  the  BI  (the  latter  are  called 
PR-Boxes).  Thus,  the  principle  of  NS  imposed  by  SR  does  not  single  out  QM  from  these  other 
post-quantum  NS  theories  [Masanes06]  (PQNS). 

These  PQNS  have  much  in  common  with  QM  such  as  no-cloning,  information-disturbance 
tradeoffs,  security  for  key  distribution,  and  others.  Recently,  van  Dam  [van  Dam05]  showed  that 
PR-Boxes  make  communication  complexity  trivial,  which  is  not  the  case  within  QM.  Other 
researchers  have  shown  that  PQNS  theories  would  lead  to  implausible  simplification  of 
distributed  computational  tasks  (see  [Pawlowski09]  and  references  therein).  It  is  now  widely 
believed  that  theories  in  which  communication/computational  complexity  is  trivial  are  very 
unlikely  to  exist.  It  is  therefore  important  to  understand  the  structure  of  the  PQNS  and  ultimately 
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to  find  physical  and  informational  principles  that  rule  them  out.  In  this  area  of  research  we  took 
steps  in  that  direction  by  investigating  the  structure  of  PR  correlations  by  forming  operators 
which  reproduce  these  PR  probability  distributions.  We  investigated  circumstances  in  which  they 
do  and  do  not  represent  valid  quantum  states. 

Bell  Inequalities  (BI) 

Nonlocality  is  expressed  by  means  of  violations  of  Bell  inequalities1  (BI)  which  set  upper  bounds 
for  classical  correlations  arising  from  local-realistic  theories.  For  bipartite  systems,  the  most  well 
know  BI  is  the  Clauser-Home-Shimony-Holt  (CHSH)  inequality7  defined  as  follows.  Consider  a 
bipartite  system  A®B,  for  Alice  and  Bob,  each  possessing  measurement  directions  A,B  =A  and 
C,D  =B  taking  measurement  values  a,b,c,d  ={±  1}.  We  define  the  correlation  E(AC)  between 
A= A  and  C=B  as 


E(AC)  =  {AC)=  £  acP(a,c  \  A,C) 

“,c={±  1}  (5) 

=  P(+,  +  140  +  P(-  -  I  A,C)~  P(+,  -  I  A,  C)  -  P{-  + 1  A,C) 


In  (5),  we  define  P(a,c\A,C)  as  the  joint  probability  that  given  the  (input)  measurement  directions 
A  for  Alice  and  C  for  Bob,  Alice  obtains  the  (output)  measurement  result  a  and  Bob  obtains  the 
value  b,  subject  to  the  normalization  condition  ’YJac={+x]P(a,c\  A,C)  =  1,  VA,C.  Finally,  we  define  the 

following  CHSH  correlation  parameter  S  by 

S  =  E(AC)  +  E(BC)  +  E(BD)-E(AD).  (6) 

S  has  been  cleverly  constructed  as  the  expectation  value  of  the  quantity  Arg  =A{C-D)  +  B(C+D). 
If  A,  B,C,D  are  classical  random  variables  taking  values  ±  1  then  it  can  be  readily  seen  that  if  (i) 
C=D,  then  |Arg|  =  \B(2C)\  =  2  and  if  (ii)  C=-D,  |Arg|  =  \A(2D)\  =  2.  Thus,  for  classical  correlation 
we  have  the  CHSH  inequality 

CHSH  inequality:  |  S  =  E(AC)  +  E(BC)  +  E(BD)  -  E(AD)  \<  Scl  =  2  (7) 

(where  the  subscript  “Cl”  denotes  “classical”).  For  a  large  class  of  measurement  directions  (but 
not  all),  quantum  states  can  violate  the  CHSH  inequality  (i.e.  |S|>2)  up  to  a  maximum  value 
shown  by  Tsirelson  [Tsirelson80]  to  be  Sq  =  2a/2.  Here,  a  quantum  state  is  defined  as  a  positive 
(i.e.  non-negative  eigenvalues)  Hermitian  matrix  with  unit  trace  denoted  by  the  symbol  p.  The 
archetypical  example  is  the  singlet  (Bell)  state 

A,*,  =(|n>-|lt»/V2  =(|01>-|1°»/V2  (8) 

with  measurement  directions  in  the  x-y  plane:  A  =  x,B  =  y,C  =  (x  +  j>)/V2  ,D  =  (x-y)/J 2  that  saturates 
the  Tsirelson  bound  with  S=  -Sq  =  -2 a/2.  This  is  a  manifestation  of  the  stronger  than  classical 
correlations  that  can  be  exhibited  by  quantum  states.  (Note:  quantum  states  with  measurement 
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directions  such  that  the  CHSH  inequality  is  satisfied,  i.e.  S< 2,  are  not  distinguishable  from 
classical  states  by  the  correlation  parameter  S). 

It  is  instructive  to  note  that  the  CHSH  inequality  in  (7)  can  be  derived  [Schumacher91]  as  a 
statement  of  a  classical  quadrilateral  inequality  for  the  correlation  metric 
A(AC)  =  i- E(AC)  =  P{+, -\A,C)  +  P(~,  +  \A,C)>0.  Substituting  this  expression  into  (7)  yields 
A(AC)+ A(BC)+ A(BD)  >  A(AD)  =>  S  =  E(AC)+ E(BC)+E(BD)—E(AD)  <  +2  (see  Figure  10).  Thus,  the 


For  a  bipartite  system  jA  ®S  with 
measurement  directions  A,B=>land  C,D=2? 

Classically:  A(AC)  +  A(BC)  +  A(BD)  >  A(AD) 

Quantum  Mechanically:  this  can  be 
violated,  e.g.  A(AC)  +  A(BC)  +  A(BD)  <  A(AD) 
with  singlet  state  (|Ti)-|tT})/.y2=(|oi)-|io))/V2 

Figure  10.  CHSH  inequality  derived  as  a  violation 

violation  of  the  CHSH  inequality  by  quantum  states  can  be  interpreted  as  a  violation  of  the 
classical  quadrilateral  inequality  which,  for  certian  measurement  directions,  yields  the  distance 
A  (AD)  via  the  direct  path  A-D  to  be  smaller  than  the  sum  of  the  distances  around  the  indirect 
path  A-C-B-D. 


'£  A 

c  A(BC)  B 


of  the  classical  quadrilateral  inequality. 


Returning  to  the  CHSH  inequality  (7),  one  notes  that  it  is  bounded  by  the  algebraic  maximum  \S\ 
<  SAm=  4.  This  follows  from  the  fact  that  the  correlations  E  are  bounded  by  |E|<1.  This  latter 
result  can  be  inferred  by  writing  E  =  P++  +  P..  -  (P+.  +  P_+)  =  2 (P++  +  P._)- 1  =  1-2 (P+.  +  P.+), 
where  =  P++  +  P__  +  P+_  +  P_+  =1  has  been  used.  Using  the  fact  that  0<  P++  +  P__  <1  and  0<  P+. 
+  P.+  <1  in  the  previous  two  expressions  for  E,  yields  the  desired  bound  |£]<1 .  Therefore,  if  the 
first  three  correlations  in  (7)  take  the  value  ±1  and  the  last  correlation  takes  the  value  +1,  we 
obtain  S=± 4.  The  implication  of  this  observation  is  that  the  regime  2x/2  <  S  <  4  represents  supra- 
correlations  that  are  stronger  than  quantum,  yet  are  unphysical  by  Tsirelson’s  bound,  i.e.  cannot 
be  realized  by  any  physical  quantum  state.  The  salient  question  to  study  is  what  ‘natural’ 
principles  determine  the  exclusion  of  such  supra-correlations.  As  a  first  hypothesis,  one  might 
surmise  that  the  principle  of  no-signaling  from  special  relativity  (i.e.  that  information  cannot  be 
instantaneously  broadcast  between  spacelike  separated  observers)  might  exclude  supra- 
correlations.  Surprisingly,  this  is  not  the  case.  In  1994,  Popescu  and  Rohrlich  (PR)  [Popescu94] 
were  able  to  construct  a  valid  joint  probability  distribution  between  a  pair  of  spacelike  separated 
observers  that  (i)  satisfies  the  non-signaling  principle,  and  (ii)  yields  the  algebraic  maximum 
correlations  allowed  by  the  CHSH  inequality.  Here  the  adjective  ‘valid’  implies  that  the  joint 
probability  distribution,  and  all  its  derived  marginal  probability  distributions,  obtain  values 
between  0  and  1,  and  satisfy  the  appropriated  normalization  requirements  (i.e.  the  joint  and  all 
marginal  probability  distributions  summed  over  all  outcomes  for  any  measurement  settings 
yields  unity).  These  correlations  are  now  called  PR  correlations,  which  we  describe  in  the  next 
section. 
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No  Signaling  (NS)  Theories  and  PR  Correlations 

We  wish  to  consider  correlations  between  n  spacelike  separated  parties  (observers)  who 

can  perform  m  possible  measurements  x\,...,x„  (x;={0,l,...,m-l),  with  r  possible  outcomes 
ai,...,an  (ai={0,l,...,r-l)  The  observed  correlations  will  be  described  by  the  joint  probability 
distribution  P{al,a2,...,an  \xl,...,xn) giving  the  probability  that  the  parties  obtain  the  measurement 
values  (outputs)  when  their  local  measurement  apparatuses  (inputs)  are  set  to  xi,. .  ,,x„. 

The  joint  probability  distribution  is  constrained  only  by  the  conditions 
0<P(al,a2,...,an\xi,...,xn)<\  and  the  normalization  condition  P{ava2,...,an  \  xp...,xj  =  1  for 

all  measurement  settings  x\ ,...ycn. 

Imposing  the  no-signaling  (NS)  constraint,  i.e.  adherence  to  the  requirement  from  special 
relativity  that  spacelike  separated  measurements  should  not  influence  each  other  due  to  the  finite 
speed  of  light  (communication),  requires  that  the  marginal  probability  distributions  satisfy  the 
additional  condition 


No  Signaling:  P(at,a2,...,ak  |Xj,...,xJ 


X  P(a1,...,an\xl,...,xn)  =  P(a1,a2 

ak+i>—>ane{  0,1} 


.,ak\xx,...,xk).{9) 


Here,  the  first  equality  in  (9)  formally  defines  the  marginal  probability  distribution  describing  the 
measurement  outcomes  of  the  first  k  parties,  when  the  last  n-k  outcomes  are  un-observed  and 
hence  summed  over.  Note,  this  marginal  probability  distribution  P{al,a1,...,ak  |x,,...,xn) formally 
depends  on  all  n  measurement  settings.  The  last  equality  in  (9)  imposes  the  NS  constraint 
requiring  that  the  marginal  probability  depends  only  upon  the  k  measurement  settings  of  the 
parties  participating  in  the  joint  measurement  (and  not  on  the  remaining  n-k  measurement  setting 
of  the  unobserved  outcomes). 

As  first  pointed  out  by  Pospescu  and  Rohrlich  [Popescu94],  the  NS  constraint  (9)  by  itself  does 
not  single  out  classical  and  quantum  theories,  i.e.  |S|  <  Sq.  PR  proposed  the  following  joint 
probability  distribution  for  two  parties  (Alice  and  Bob)  with  two  measurement  settings  (inputs) 
x,y  ={0,1},  and  two  measurement  outcomes  (outputs)  a,b  ={0,1}  given  by 


PR  Box:  P(a,  b  \  x,  y)  = 


if  a  ®  b  =  x-y 
otherwise 


(10) 


By  considering  all  possible  inputs  and  outputs,  it  is  straightforward  to  show  that  PR  correlations 
of  (10)  satisfy  all  the  requirements  for  a  NS  theory  as  follows:  normalization  (total  probability) 
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Y  P(a,b\x,y) 

a,b={  0,1} 

=  P(0,0|x,>’)  +  P(l,l|x,>’)  +P(0, 1 1  x,  y)  +  P(l,  0 1  x,  7) , 
a  ©Z?  =  0  a  ©Z?  =  1 

(1/2  +  1/2)  (1/2  +  1/2)^, 

~  ^0  ,x-y  ^1  ,x-  v’  ’ 

=  1  Vx,y, 


and  the  NS  constraint 


P(a\x,y)=  Y  P(a,b\x,y ) 

*={0,1} 

=  P(fl,0|  x,y)  +  P(g,l|x,y) 

a®b=a®0=a  a®b=a®\=a 

=  \l2SatX.y+\l2  SS>x.y,  (12) 

Jl/2  +  0  (if  a  =  0  &  x  •  j;  =  0),  0  +  1/  2 (if a  =  0&x-y  =  l) 

{0  +  l/2(ifa  =  l&x-y  =  0\l/2  +  0(ifa  =  l&x-y  =  l)  ’ 

=  1/2  \/a,x,y, 

=  P(a  |  x)  \/a,  x  ( =>  Isotropic ,  i.e.  P(a  |  x)  =  1  /  2  indep  of  a,  x). 

With  the  PR  Box  define  above  in  (10)  we  can  compute  correlations  as 

E(a,  b\x,y)  =  P(0, 0 1  x,  y)  +  P(  1, 1 1  x,  -^(0, 1 1  x,  y)  -  P(l,  0  |  x,  y) , 

a®b  =  0  a®b  =  1 

(1/2  +  1/2)  b^Xy  - (1  / 2  +  1/2)^,  (13) 

=  |+1  if  x-j  =  0,  i.e.(x,y)  e  {(0,0),  (0,1),  (1,0)}, 

[-1  if  x-y  =  1,  i.e.(x,y)  =  (1,1), 

where  we  have  used  E(a,b  |  x,y)  =  ^a,6,e{+ l}a'MP(a,b  \  x,y),  where  a'  =  1-2 a  (//  =  1-2Z>)  associates  the 

measurement  values  a'(b')  e  {+1,-1}  with  the  measurement  value  labels  (bits)  a(b)em\, 
respectively.  Therefore,  in  Figure  10,  assigning  Alice’s  measurement  directions  A,B  =A  the  bit 
labels  x^=l  and  xg=0,  and  Bob’s  measurement  directions  C,D  =B  the  bit  labels  vc=0  and  jy>=  1, 
and  using  (6)  yields  the  algebraic  maximum  Sm=A  of  the  CHSH  inequality,  as  illustrated  in 
Figure  11. 
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PR  Box: 

P(a,b\x,y)  = 


if  a®b  =  x-y 
otherwise 


with  measurements  x,y  and  outcomes  a.b  as  bits, 
i.e  a,b,x,y  e  {0,1}  (Note:  {0,1}  <-»  {+1,-1}) 


(vc=0)C  A(BC)  BCYs=0) 


PR  Box  yield  algebraic  maximum  of  S 

5  =  E(AC)  +  E(BC)  +  E(BD)  -  E(AD) 
=  1  +  1  +  1  -  (-1) 
=  4 


Figure  11.  PR  Box  with  joint  probability  distribution  achieving  the  algebraic  maximum 
Sm= 4  of  the  CHSH  inequality. 


Since  Sm=4>Sq=2^2,  no  quantum  (i.e.  physically  realizable)  state  can  reproduce  the  above  PR 
probability  (10).  However,  the  following  “state”  [AcinlO]  0  =  ^+|o+^0+|  +  ^_|o_^0“| with 

a±  =(l±V2)/2and  Bell  states  |o±)  =  (|00)  +  |ll)>/V2,  yields  the  PR  probability  (10)  through  the  usual 
trace  rule  PPR(a,b\x,y)  =  Tr[OMxa  ®My]  with  =  and  {M?,MxbD}  =  {(O',  +<T2)/V2, 

(<7,-<j2)/V 2},  where  {o-,}ie{123}are  the  usual  Pauli  matrices.  Note  that  the  form  of  the  joint 
measurement  between  Alice  and  Bob  written  as  a  pure  tensor  product  of  local  observables 
Mxa  ®Myb,  ensures  the  locality  of  the  spacelike  separated  measurements,  which  cannot  increase 
entanglement  between  the  parties.  A  measurement  involving  the  sum  of  pure  tensor  products, 
such  as  Mx  ®  Ml  +  M'x  ®  M’-r  which  might  possibly  create  entanglement,  would  involve  non-local 
measurements  between  the  parties,  which  could  only  be  physically  realized  if  the  parties  were 
brought  together.  The  important  point  is  that  O  does  not  represent  a  physical  quantum  state  since 
it  is  non-positive,  i.e.  it  possesses  the  negative  eigenvalue  a_  =(l-V2)/2.  Henceforth,  we  shall 
refer  to  non-positive,  unit  trace  Hermitian  operators  O  capable  of  producing  NS  probability 
distributions  as  “states,”  and  reserve  the  specific  term  “quantum  state”  or  “q-state”  for  the 
physically  realizable  positive,  unit  trace  Hermitian  operators  denoted  as  p>0,  (i.e.  density 
matrix). 


Following  Acin  et  al.  [AcinlO]  we  desired  to  investigate  all  sets  of  n-party  spacelike  correlations 
in  terms  of  local  quantum  observables  (measurements)  Mnon_sig  =M*>  ®  •  •  •  ®  Mx,‘  that  ensure  NS. 
These  correlations  can  be  written  in  the  form 


P0=Po(aV-’an\XV-’Xn)  =  Tr[°Ma' 


(14) 


Without  loss  of  generality,  we  can  take  the  local  measurement  operators  Mxa  =  IT  =|a) ■  (a|to  be 

the  projection  operators  onto  “spin-component”  a  in  the  “direction”  x.  Requiring  that  proper 
probabilities  be  derived  from  all  local  quantum  measurements  imposes  the  condition  that  O  be 
positive  on  all  product  states.  This  implies  that  0=W  is  an  entanglement  witness  (EW,  see 
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[Guhne09])  with  the  property  ( a,p,...\w\a,p,..)>0 .  Here  some  definition  are  helpful.  A  q-state 
is  separable  (contains  only  classical  correlations)  if  it  is  of  the  form  psep  =  'YJjpip?'  ®P-2  ®  •  •  •  ®  pfN 
where  each  pf'  is  a  local  density  matrix  and  ^  p,  =  l.  (If  a  q-state  is  not  separable,  it  is  entangled). 
Each  local  density  matrix  has  a  (non-unique)  ensemble  decomposition  Pi*  =ZjPlH)H\  where 
^jPij  =1.  The  requirement  that  W  is  positive  on  all  product  states  (a,p,...\ W | a, /?,...)>  0 ensures 

that  Tr\p  'vw\  >  o  from  the  form  of  psep.  A  q-state  p  such  that  Tr\pw\  <  o  is  then  entangled  (since  it 
is  not  separable),  and  W  is  said  to  “witness”  (or  exhibit)  the  entanglement  of  p.  Note  that  W  is  in 
general  a  non-positive  Hermitian  operator.  In  the  context  of  (10),  we  now  consider  O— as  a 
state  (not  necessarily  a  q-state)  from  which  to  derive  NS  correlations  through  the  joint 
probability  distributions 

Pw  =  P(ax,...,aN  \  x\,...,xn)  =  Tr[W M*1  ®---®Mxa"  ]  >  0.  (15) 

The  correlations  (4)  are  termed  Gleason  correlations  by  Acin  et  al.  [AcinlO], 

The  subtle  distinction  between  (14)  and  (15)  is  that  the  latter  produces  positive  probabilities  for 
all  local  NS  measurements,  while  the  former  may  produce  non-negative  probabilities  on  only  a 
subset  of  NS  measurements.  This  distinction  is  important  since  it  has  been  shown  [Guhne09, 
BamumlO]  that  for  bipartite  systems  n= 2,  any  Gleason  correlation 
P(al,a2 1  x,,x2)  =  Tr[W M*1  ®M*2 ]  >  0  can  be  converted  to  a  probability  distribution  derived  from  a  q- 

state  p^fb)  =  |oM)(ops |  with  modified  measurements P(at , a2  | xi,x2)  =  Tr[WM[ '■  ® ] 
=  7>f/y,Vs .  M;'1  ®  ]  >  o.  Here  |Ora)is  any  pure  bipartite  state  (not  necessarily  maximally 

entangled).  The  proof  relies  on  the  explicit  use  of  the  Choi-Jamiolkowski  isomorphism  (CJI) 
[Guhne09,  BamumlO,  Vedral97]  which  allows  any  bipartite  (n=2)  witness  W  to  be  written  as 
W,n=2>  =(/ ®  A)(p0n}),  where  A  is  a  positive  trace  preserving  map.  In  the  above,  M*2  =A *(M£) 

where  A*  is  the  adjoint  of  the  map  A,  i.e.  Tr[AA(5)]=Tr[A*(^)5],  The  proof  then  follows  directly  as 
Pw(a,b  \x,y)  =  Tr[W  Mxa  ®Myb}  =  Tr[{I®  A)(p,  ,)Ma*  ®M£] 

—  (16) 

=  Tr[M*a  ®  Myb  (7  0  A)(/7|0^))]  =  Tr[Mxa  ®  A*(A/')^*„)]  =  Tr[p^p)  M*a  0M/)], 

where  the  second  equality  uses  the  CJI,  the  third  equality  uses  the  cyclic  property  of  the  trace, 
the  fourth  inequality  utilizes/®  A  acting  to  the  left  on  the  tensor  product  measurements  ,v/;  ®Ml 
thereby  introducing  the  adjoint  A*  operation  and  the  modified  local  measurement  operation 
=A*(M*2)in  the  last  equality.  Acin  et  al.  [AcinlO]  point  out  that  the  CJI  decomposition 

Wm2)  =  (/®A)(pjiVi  /)  in  general  fails  for  n>2  (which  they  demonstrate  by  a  specific  example). 

Thus,  the  Gleason  correlations  (4)  are  strictly  larger  (|A|  >Sq)  than  quantum  correlations  for  n>2 
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(and  equivalent  only  for  n< 2).  The  state  <9  =  a+|<t>+^®+|  +  a_|<t>  ^<t>  |used  in  the  example  of  PR 

correlations  in  the  discussion  after  Figure  11  is  not  an  EW  since  it  can  produce  negative 
probabilities  for  measurements  other  than  those  considered  (it  would  be  an  EW  if  it  produced 
positive  probabilities  for  all  measurement  choices).  Acin  et  al.  [Acini 0]  classify  the  distributions 

P(a],...,an  |  x, . xn ) as  (i)  No-Signaling  if  and  only  if  P  can  be  written  in  the  form  of  (14),  (ii) 

Quantum  whenever  O  is  positive  (O>0),  and  (iii)  Local  if  and  only  if  O  corresponds  to  a 
separable  quantum  state.  In  the  following,  we  investigate  the  NS  correlations  of  (10)  and  the 
conditions  for  which  they  become  either  Gleason,  or  Quantum  correlations. 

4.0  RESULTS  AND  DISCUSSION 

4.1  Multipli-entangled  photons  from  a  spontaneous  parametric  down-conversion  source 

Experimental  analysis  and  testing  apparatuses  for  Schioedtei  are  very  similar  to  that  for  any 
SPDC  source.  With  the  more  complex  ring  pattern  generated  though  there  are  modifications  one 
must  do  to  the  standard  detection  scheme.  The  experimental  configuration  for  Schioedtei  is 
shown  in  Figure  12.  The  testbed  consists  of  a  violet  (405  nm)  femtosecond  pulsed  pump  source 
(Millenia  PRO  15sJ  >  Tsunami  3960- 15HP  >  Inspire  Blue  FM)  with  an  average  power  of  ~1.4 
W,  ~100  femtosecond  pulses  and  a  repetition  rate  of  80  MHz.  The  405  nm  pulses  first  pass 
through  a  ~  12.5  mm  quartz  pre-compensator  and  a  half-wave  plate  set  to  22.5°  to  rotate  the  input 
linear  polarization  to  the  required  45°  for  equal  excitation  of  the  crystals  before  entering 
Schioedtei.  Proper  alignment  of  the  crystal  was  accomplished  with  live  images  from  a  cooled 
CCD  camera  (Princeton  Instruments  Pixis  1024BR).  The  photons  were  collected  in  free  space 
collimators  located  1.5  meters  behind  Schioedtei.  This  distance  is  the  minimum  amount  required 
to  obtain  the  useable  spatial  separation  required  for  detector  access  to  the  middle  blue  diamond 
of  intersection  points  (5,  6,  7,  and  8).  The  post-compensating  crystals,  inserted  in  the  down- 
converted  photon  paths,  are  8x8x1  mm  type  II  phase  matched  P-BBO  (at  angles  of  0  =  41.9°  and 
cp  =  30°)  as  Schioedtei ’s  orientation  is  non-collinear  and  there  is  no  interaction  between  the  pump 
and  the  compensators.  These  compensators  could  not  be  used  for  compensation  of  a  collinear 
configuration  as  they  were  phase  matched  for  SPDC  at  810  nm  when  exposed  to  a  405  nm 
excitation  beam.  Photon  collection  was  accomplished  via  fiber  coupled  collimators  immediately 
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Figure  12.  Experimental  testbed  to  analyze  the  Schioedtei  source. 

followed  by  2  nm  bandpass  filters.  The  output  of  the  bandpass  filter  was  routed  directly  into 
fiber-coupled  single  photon  counting  avalanche  photodiodes  (APDs)  (Perkin  Elmer  SPCM- 
AQ4C).  Coincidence  detection  was  accomplished  by  a  four  channel  coincidence  counting 
module  (CCM)  [Branningll], 


Figure  13.  False  color  CCD  images  of  custom  crystal  assembly  (1  sec  exposure).  A,B  are  the  type  II  non- 
collinear  outputs  from  each  individual  crystal.  C  is  the  combined  output  from  the  crystal  stack. 

A  trio  of  false  color  CCD  camera  images  of  Schioedtei  output  is  shown  in  Figure  13.  The  twelve 
overlap  regions  are  clearly  visible  and  the  spatial  symmetry  of  the  output  should  be  clearly  noted. 
The  orientation  of  the  crystal  assembly  gives  an  approximate  Gaussian  profile  on  spots  5, 6, 7, 8 
and  a  slightly  elongated  profile  for  spots  1,2,3,4,9,10,11,12.  The  alignment  image  in  Figure  13  is 
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utilized  for  aligning  the  proper  orientation  of  the  rings  while  a  back  propagated  beam  shown  in 
Figure  14  aligns  the  collimators  to  the  intersection  points. 


Figure  14.  Alignment  image  of  the  Schioedtei  crystal  stack. 

As  stated,  Schioedtei  was  constructed  from  P-BBO  though  any  type  II  material  can  be  used. 
Materials  such  as  BiBO  (Bismuth  Borate,  BiB306)  have  been  shown  to  have  a  higher  photon 
generation  rate  than  P-BBO  [Rangarajan09]  and  this  will  be  the  next  step  for  Schioedtei. 
Secondly,  increasing  the  useable  photon  count  rate  in  Schioedtei  can  be  accomplished  by 
factoring  the  GVM  phase  matching  constraint  [U’Ren06]  into  the  crystal  construction.  A  GVM- 
matched  configuration  [FantolO]  is  possible  by  alternating  reduced  thickness  Schioedtei  and  a- 
BBO  layers.  a-BBO  can  be  used  as  a  compensator  since  there  is  no  second  order  nonlinear  effect 
in  a-BBO  crystal  due  to  the  centric  symmetry  in  its  crystal  structure.  Such  a  GVM  source  would 
provide  the  same  up  to  six  spatially  separate  entangled  pairs  as  Schioedtei,  while  alleviating  the 
need  for  spectral  filtering  of  the  photons.  An  increase  in  useable  signal  rates  of  lOx  over  a  typical 
type  II  source  is  realizable  with  GVM  matching. 

Schioedtei  source  uses  and  applications 

Another  applicable  area  of  extreme  interest  is  in  the  generation  of  photon-based  cluster  states. 
Cluster  states  play  a  central  role  in  the  measurement-based  one-way  quantum  computation 
approach  [Walther05,  RaussendorfOl].  In  this  scheme,  the  entanglement  resource  is  provided  in 
advance  through  an  initial,  highly  entangled  multi-qubit  cluster  state  and  is  consumed  during  the 
quantum  computation  by  means  of  single-particle  projective  measurements.  The  feedforward 
nature  of  the  one-way  computation  scheme  renders  the  quantum  computation  deterministic,  and 
removes  much  of  the  massive  overhead  that  arises  from  the  error  encoding  used  in  the  standard 
quantum  circuit  computation  model  [O’Brien07].  Figure  15  illustrates  a  scheme  for  utilizing  the 
output  of  Schioedtei  to  generate  a  four  photon  cluster  state,  |C)4  [Schmid07].  This  particular 
example  employs  the  spots  1,2, 3,4  and  requires  insertion  of  two  half-wave  plates,  a  SWAP  gate 
and  a  controlled-phase  (CPhase)  gate.  This  scheme  could  be  expanded  to  include  the  other  eight 
spots  to  generate  even  larger  cluster  states.  Such  experiments  are  currently  being  explored  in- 
house. 

More  complex  cluster  states  can  be  constructed  from  Schioedtei  with  additional  hardware.  This 
includes,  but  is  not  limited  to,  the  construction  of  box  cluster  states  [Prevedel07,  Walther05],  In 
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fact  Schioedtei  is  capable  of  producing  two  4-qubit  box  states  simultaneously  by  using  eight  of 
the  spots:  1,2, 3, 4  and  5, 6, 7, 8.  As  the  states  Schioedtei  outputs  at  these  two  sets  of  spots  are 
different, 


f 

V 


Hadamard 
Compensator 
X  Rotation 


SWAP  Gate 


C  Phase  Gate 


I c)4  =  \  (\HHHH)wa  +  \HHW)wa  +  \VVHH)wa  -  \VVVV)WA) 

Figure  15.  Experimental  setup  for  4-qubit  cluster  state  generation  utilizing  Schioedtei. 

slightly  different  preparation  methods  are  required  for  the  two  boxes,  as  shown  in  Figure  16. 
After  the  preparation  is  complete  the  two  box  states  are  completely  equivalent.  With  additional 
preparation  and  resource  photons  these  states  can  be  used  as  the  building  blocks  of  larger  states 
such  as  the  6-qubit  butterfly  network  [MalO,  SoedalO]. 


1 

3 
2 

4 


|^)bOX—  Swap2,3  H1H2H3H4  CZ2,4  Swap2,4  X3X4|l/;)i)2,3,4 
|*/>)bOX=  Swap6,7  H5H6H7H8  026,81^)5,6,7,8 

Figure  16.  Experimental  construction  of  a  4-qubit  box  cluster  state  utilizing  Schioedtei. 

An  advantage  of  the  Schioedtei  configuration  is  the  diversity  of  states  that  it  is  capable  of 
generating.  Schioedtei  allows  for  the  direct  generation  of  the  (unnormalized)  state  |  HV)  ± 
el(p\  VH)  along  with  the  generation  of  the  state  \HH)  ±  el(p\VV)  with  the  addition  of  a  half-wave 
plate.  In  addition,  separable  states  such  as  | HV)  ±  el<p\VV)  or  | HV)  ±  el(p\HH)  can  also  be 
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directly  generated  with  clever  combinations  of  the  twelve  output  intersections  and  proper 
compensation. 

4.2  A  multi-layer  three  dimensional  superconducting  nanowire  photon  detector 
General  amplitude  amplification 

We  now  take  a  closer  look  at  the  minimum  three  layers  needed  to  create  the  device,  shown  in 
Figure  17.  The  bottom  layer  17a,  consists  of  non-superconducting  leads  (red)placed  on  an 


Figure  17  A  plain  view  of  the  three  layers  in  the  multilayer  design. 

insulating  substrate  (gray),  such  as  R-plane  sapphire,  MgO  or  Si.  Note  that  there  is  no  complete 
circuit  on  this  level,  so  the  current  will  be  forced  to  move  up  to  the  next  level.  The  optimal 
minimal  spacing  will  depend  on  the  insulating  ability  of  the  substrate  to  prevent  leakage  and 
crosstalk,  mainly  between  the  input  and  output  channels  but  also  with  the  superconducting 
nanowires  passing  above.  Over  the  bottom  layer  will  be  a  second  layer  of  deposited  substrate 
17b.  This  layer  will  then  have  vias,  i.e.  holes  (green),  which  pass  completely  through  it  (gray)  at 
predetermined  locations  so  as  to  hit  the  input/output  leads  in  the  bottom  layer.  These  vias  are 
fdled  with  superconducting  material  (in  practice  it  may  be  advantageous  to  use  non¬ 
superconducting  material  here,  if  the  pixels  are  long  enough  to  avoid  the  latching  condition)  thus 
completing  the  middle  layer  of  the  device.  Alignment  will  be  a  very  important,  but  not 
insurmountable  issue,  as  these  structures  are  on  the  scale  of  approximately  100  nm  in  width  and 
current  alignment  techniques  can  achieve  results  on  the  order  of  1  nm  [Anderson04],  Finally  the 
detection  layer  Figure  17c,  will  be  deposited  on  top  of  the  middle  layer.  Alignment  of  the 
superconducting  bridges  with  the  vertical  “posts”  in  the  middle  layer  will  be  important  for  the 
overall  detection  efficiency  [Kerman07]. 
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The  device  will  have  significantly  higher  number  resolution,  while  maintaining  a  useful 
detection  area.  It  has  several  parameters  which  can  control  the  reset  time  to  avoid  latching  while 
still  minimizing  the  rest  time.  An  array  of  pixels  of  arbitrary  number,  size  and  shape  is  possible. 
Most  of  the  detector  will  remain  active  after  a  single  photon  is  absorbed  as  opposed  to  small 
number  or  single  meander  detectors  which  are  effectively  blinded  by  a  single  photon.  The  active 
area  of  the  detector  can  be  tuned  by  changing  the  number  or  the  shapes  of  the  pixels.  These 
advantages  are  compelling  theoretical  evidence  for  the  construction  and  testing  of  multi-layer 
superconducting  number-resolving  photon  detectors. 


Figure  18  shows  a  toy  model  of  a  working  multi-layer  SNSPD.  For  clarity  only  the  super-cooled 
part  of  the  chip  is  shown.  The  dark  green  substrate  and  red  leads  are  similar  to  the  grey  substrate 
and  red  leads  in  Figure  17a.  The  blue  vertical  posts  are  shown  in  a  semi-transparent  middle 
layer  similar  to  Figure  17b  for  clarity.  The  blue  superconducting  pixels  on  the  surface  are 
similar  to  Figure  17c.  The  development  of  this  design  has  resulted  in  two  patents  filled  by  the 
AFRL  JAG  officer  with  the  U.S.  patent  office. 

4.3  Laboratory  upgrade  and  ongoing  research  in  integrated  waveguide  quantum  circuits 

The  Quantum  Information  Science  Laboratory  originally  located  in  lab  18  in  Bldg.  104  was 
relocated  to  a  larger  facility  in  Bldg.  3  Suite  15.  The  transition  to  the  new  facility  allowed  for  the 
addition  of  2  more  optical  tables,  multiple  work  benches  and  equipment  storage  cabinets.  The 
facility  is  partitioned  in  two  separate  work  areas  defined  by  laser  curtains  allowing  separate 
experiments  to  be  conducted  concurrently  with  lasers  of  class  4  or  lower.  During  this  period  the 
Ti:Sapphire  laser  was  upgraded  from  a  femtosecond  1.5  W  system  to  a  femto/picosecond  3.5-4 
W  system.  Additionally  installed  were  both  femtosecond  and  picosecond  second  harmonic 
generation  units  (SHG)  were  added  to  the  system  to  generate  powers  greater  than  1  W  in  the 
blue/violet  regime.  These  additions  completed  the  upgrade  to  the  entangled  photon  generation 
testbed. 
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Further  effort  has  been  placed  to  reduce  the  footprint  size  of  quantum  gates/circuits  built  from 
bulk  optical  components.  This  added  research  focuses  on  the  use  of  integrated  optical 
waveguides  to  construct  the  quantum  gates/circuits.  The  direction  of  the  research  exploits  two 
arenas:  (i)  world  class  domestic  researchers  at  Rome  Research  Site  and  WPAFB  along  with 
universities  such  as  Columbia,  MIT  and  RIT,  and  (ii)  and  world  class  international  researchers 
through  EOARD  at  universities  such  as  Bristol,  Oxford  and  Vienna.  Expanding  the  ongoing 
research  in  optical  waveguides  was  a  necessary  step  and  made  possible  in-house  with  the 
acquisition  of  an  optical  wafer  probe  station.  The  probe  station  along  with  multiple  table  top 
probe  stations  will  be  utilized  for  the  testing  and  integration  of  quantum  photonic  integrated 
circuits  (QPIC).  The  entangled  photons  generated  by  the  existing  generation  testbed  are  routed 
into  the  QPICs  to  validate  the  chips  functionality.  The  acquisition  of  a  second  Ti:  Sapphire  laser 
and  optical  parametric  oscillator  (OPO)  expanded  the  testbed’s  available  wavelength  range  from 
the  original  span  of  600-1000  nm  to  a  span  of  340-2500  nm.  The  OPO  greatly  increases  the 
diversity  of  materials  that  the  QPICs  can  be  constructed  from.  The  upgraded  components  have 
arrived  and  the  full  testbed  is  under  construction. 

4.4  Theory/experimental  requirements  of  imperfect  two-qubit  linear  optical  photonic  gates 

The  optimization  method  we  have  developed  maximizes  the  success  probability  S  for  a  given 
target  transformation  ATar,  for  given  ancilla  resources,  and  for  a  given  fidelity  level  F<1.  This  is 
mathematically  equivalent  to  unconstrained  maximization  of  the  function  S+F/e  in  the  space  of 
all  matrices  U,  where  1/e  is  a  Lagrange  multiplier.  Here  e — >0+  corresponds  to  maximizing  the 
success  probability  while  requiring  perfect  fidelity  (F  =  1).  As  e  is  increased,  the  maximum  of  S 
+  F/e  yields  linear  optics  transformations  that  maximize  the  success  S  as  a  function  of  the  fidelity 
F.  Given  one  transformation  U  that  (locally  or  globally)  maximizes  success  S  for  a  given  fidelity 
F,  e  may  be  continuously  varied  to  obtain  a  one-parameter  family  of  optimal  transformations, 
tracing  out  a  curve  in  success-fidelity  space.  Note  that  in  general  the  members  of  these  families 
need  not  be  all  unitary,  however  for  some  gates  of  interest,  including  the  CZ  gate,  all  members  of 
the  family  are  unitary.  Figure  20  shows  optimal  results  for  the  CZ  gate.  Here  each  point 
corresponds  to  a  unique  unitary  mode  transformation  U.  As  previously  reported  we  find  an 
interesting  feature  of  these  unitary  matrices.  The  optimal  solution  with  fidelity  F  =  1  was  found 
by  Knill  to  have  a  surprising  form  [Knill02],  which  we  have  dubbed  the  “Knill  Form" 
[Uskov09],  where  one  mode  of  each  qubit  is  non-interacting,  e.g.,  in  the  CZ  case  U  acts  as  the 
identity  on  modes  1  and  3  (or  equivalently  1&4,  2&3,  or  2&4).  This  form  has  been  found  to  hold 
for  the  CZ  gate  and  for  the  TS  Toffoli  Sign  gate  (CNOT  and  Toffoli  respectively  are  equivalent 
to  these  up  to  local  rotations). 

We  now  propose  an  experiment  that  will  test  the  results  shown  in  Figure  19.  Reck  et  al.  have 
shown  that  any  discrete  N  x  N  unitary  transformation  U  can  be  implemented  as  a  multi-port 
device  consisting  only  of  variable  transmittance  beamsplitters  and  phase  shifters  [Reck94],  Their 
method  is  a  decomposition  in  which  each  unitary  matrix  element  below  the  diagonal  is 
transformed  into  zero  by  a  2  x  2  rotation  matrix  embedded  in  an  N  x  N  matrix  which  is  otherwise 
equal  to  the  identity.  For  example,  the  2  x  2  rotation  acting  on  modes  N  and  N-l,  which 
eliminates  the  element  Un,n-i,  takes  the  form  Tn  n_i shown  in  Figure  19. 
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Tn,n- i  — 


0  \ 

el<^sin(a>)  e^sin^a))  I 
e^sin (o>)  e^sin(co)/ 


Figure  19.  Rotation  matrix  for  modes  N  -1  and  N. 

The  method  is  recursive  and  requires  one  iteration  for  each  pair  of  modes.  Finally,  we  obtain 
U(N)Tnn_1Tnn_2  —Tz.iD  =  I  where  D  is  a  diagonal  matrix  of  phases.  The  desired 
transformation  U  is  then  decomposable  as  U(N)  =  D-1'r2“11T3“11  ...  T^_±.  Physically,  each  NxN 
transformation  7)”1  is  implemented  as  a  variable  transmittance  beamsplitter  with  a  phase  plate  on 
one  input  mode,  while  D1  corresponds  physically  to  a  phase  shift  on  each  output  mode 
[Reck94],  Thus  a  generic  two-qubit  operation,  which  needs  at  least  N  =  7  modes  (Nc  =  4 
computational  modes  and  Na  =  3  ancillas)  requires  a  minimum  of  21  beamsplitters  and  28  phase 
shifters.  A  controlled  unitary  gate  (N  =  Nc+Na  =  4+2  =  6)  requires  at  least  15  beamsplitters  and 
21  phase  shifters.  If  unitary  dilation  is  required  (as  is  often  the  case)  the  number  of  optical 
elements  increases  rapidly.  However  our  experiment  does  not  require  unitary  dilation  and 
furthermore  as  noted  by  Reck  et  al.,  if  an  element  of  the  unitary  matrix  is  already  zero,  then  no 
transformation  is  required.  The  element  is  skipped. 

Success  Rate  of  the  CNOT  gate  vs.  8 


Figure  20.  Improved  success  rates  for  compromised  S. 


Here  we  return  to  the  “Knill  Form,"  where  in  the  case  of  CZ  we  find  that  nine  of  the  elements 
below  the  diagonal  are  already  zero.  Therefore  the  unitary  transform  can  be  implemented  with 
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only  six  beamsplitters  and  ten  phase  shifters.  We  can  perform  this  decomposition  for  each  data 
point  in  Figure  20,  and  find  the  rotation  angles  coy  and  phases  4>ij  in  each  case.  Surprisingly  we 
find  numerically  that  all  of  the  phase  shifts,  <J>ij ,  are  constant  along  the  entire  length  of  the  curve 
in  Figure  20.  Therefore  only  the  six  beamsplitter  rotation  angles  coy  out  of  a  total  of  36  possible 
variables  need  to  be  modified  to  vary  ,  making  the  experiment  much  more  physically  realizable. 
To  be  specific,  the  transformation  only  requires  beamsplitters  acting  on  the  following  mode 
pairs:  (i;  j)  =  (6;  5);  (6;  4);  (6;  2);  (5;  4);  (5;  2);  (4;  2).  Figure  21  shows  that  the  six  beamsplitter 
rotation  angles  change  smoothly  with  6.  Implementing  such  rotations  and  constant  phase  shifters 
will  recreate  the  unitary  matrices  from  Figure  20. 


1-F 


Figure  21.  Beamspliter  transitivities 

This  system  lends  itself  to  being  implemented  with  2x2  Mach-Zehnder  interferometers  (MZI)  in 
place  of  standard  beamsplitters.  The  transmittance  of  the  MZI  is  controlled  dynamically  by 
adjusting  the  phase  difference,  without  having  to  alter  the  physical  system.  These  interferometers 
have  already  been  put  on  optical  chips  by  Thompson  et  al.  [Sohma94]  among  others.  Indeed, 
significantly  larger  electro-optical  matrix  switches  have  been  proposed  and  built  for  broadband 
optical  communication  networks  [Sohma94,  Drever83],  Figure  22  shows  a  multi-port  device  that 
mixes  seven  input/output  modes  (thin  lines)  using  2x2  variable  transmittance  beamsplitters 
(rectangles),  each  of  which  has  a  phase  shifter  on  one  of  its  input  modes  (ellipses).  An  additional 
phase  shifter  is  placed  on  each  device  output  mode.  The  thick  line  is  a  simple  mirror.  J.  L. 
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O'Brien  recently  proposed  a  similar  7x7  single-chip  MZI-based  device  made  from  lithium 
niobate  waveguides  [Sohma94],  The  intended  purpose  of  this  chip  was  to  be  able  to  perform  any 
two-qubit  unitary  operation,  i.e.  any  transformation  in  SU(4).  However,  such  a  device  would  also 
be  capable  of  performing  the  experiment  described  above. 


Variable  Tran snrwttance  Beamsplitter  -  Optical  Made 

Phase  Shifter  Mirror 


Figure  22.  General  multiport  device  schematic. 

4.5  Nonlocality,  entanglement  witnesses  and  supra-correlations 
No  Signaling  (NS)  Correlations:  2-Qubits 

Following  Acin  et  al.  [AcinlO]  we  define  an  » -partite  probability  distribution a _  | 
as  being  NS  if  and  only  if  there  exists  local  quantum  measurements  and  a  Hermitian  operator 

O  of  unit  trace  such  that  (3)  holds.  It  is  important  to  note  that  O  need  not  produce  positive 
probabilities  for  other  measurements  outside  this  set.  Acin  et  al.  [AcinlO]  give  a  prescription  for 
the  formal  construction  of  O  given  the  set  of  measurements  M*1 .  In  the  following  we  present  an 

explicit  construction  for  O  for  the  case  of  n= 2  qubits  (r= 2  outputs,  i.e.  a,b  =  {0,1})  and  arbitrary 
number  m  of  measurement  inputs  (x,y  =  {0,1,..., m-1 }).  Later,  we  extend  this  to  the  case  of  n= 3 
for  qubits. 

As  stated  in  Section  3.4,  without  loss  of  generality  we  can  take  the  local  Hermitian  measurement 
operators  to  be  the  projection  operators  onto  “spin-component”  a  in  the  “direction”  x, 
Mx  =nxa  =|a)  (a|.  For  each  x,  the  completeness  of  the  measurement  operators  give 
where  lrxr  =  1  is  the  rxr  identity  matrix.  This  allows  us  to  write  the  a=r-l  measurement  operator 
as  Mxa=r_1=lrxr-YJ'~^0Mxa.  One  defines  the  (tilde)  Hermitian  matrices  Mx  dual  to  Mx through  the 
inner  product  Tr[MxMxJ]  =  Sxx,Saa..  For  the  bipartite  case  n— 2,  with  in  general  m  measurement 
settings  with  r  measurement  outcomes,  one  has 


r-2  m- 1 


r-2  m- 1 


r-2  m- 1 


o  =  Y,  Z  p(a’b  l  x’y)Mxa  |  x)Mxa  0  7  +  X  Y, PQ>  I  y)i  ®  Myb  + h ®  7,  ( : 1  ■ 7) 


a,b= 0  x,y= 0 


a= 0  x=0 


b= 0  y 
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where  7 is  the  tilde  matrix  dual  to  the  rxr  identity  matrix  /,  with  the  additional  orthogonality 
conditions  defined  by  7>[/  7]  =  7>[  7]  =  l,  Tr[MxJ]  =  0,  and  Tr[iMxa',]  =  Tr[Mxa',]  =  0.  The  conditions 


Alice  and  Bob  share 
paired  (generalized) 
PR  Boxes  with 

P(a,b\x,y)  =  \/2 
iff  a@b  =  x-y  mod 2, 
(0  otherwise) 


P(a,b\  x,  y)  =  Tr[0 Mxa  ®Mjv]>0 


Figure  23.  PR  Box  shared  between  Alice  and  Bob. 


ensure  that  O  is  Hermitian,  Tr[0]=l  and  probabilities  are  given  by  the  trace  formulas 
P(a,b  \x,y)  =  Tr[OM *  0 Myb  ] ,  P(a  \  x)  =  Tr[0 Mxa  ®  I]  and  P(b  \  y)  =  Tr[Oi®Myb}.  This  is  illustrated  in 
Figure  23  where  Alice  and  Bob  share  PR  correlations  by  means  of,  what  are  termed  in  the 
literature,  a  pair  of  PR  boxes  (or  NS  {non-signaling}  boxes). 

In  the  following  we  specialize  to  the  case  of  qubits  (r= 2,  a,b=  {0,1 })  with  arbitrary  number  in  of 
measurement  inputs  (x,y  =  {0,1,..., m-1  }) .  In  this  case  the  measurement  operators  m;=ii  are  given 
as  projection  operators  for  “spin-up”  along  the  directions  x^mx  on  the  Bloch  sphere.  The  m;=ii 
are  just  density  matrices  on  the  Bloch  sphere  written  as 

MLo  =  \  °),  (0 1  =  l/2  i1  +  mx  ■  d-),  mx  =  mx  (sin  0x  cos  <j>x ,  sin  0X  cos  0X ,  cos  0X  ) 

\mx  \  <  1,  (density  matrix  on  Bloch  Sphere), 

where  <r  =  (cr1,cr2,cr3)is  the  vector  of  single  qubit  Pauli  matrices.  Although  not  required  for  the 
case  of  qubits,  the  projection  onto  “spin-down”  along  x  is  given  by 
=14(11  =  1/2(7  - mx  ■a)  =  I — Mxa=0 ,  with  /  the  2x2  identity  matrix.  Equation  (17)  now  simplifies 
to  the  form 


m— 1  m— 1  m— 1 

0=  X  P(a  =  0,b  =  0\x,y)M’®M’+Y,P(a  =  o|  x)m;®i  +  Zp<b  =  0|y)7  0Mq  +7  07.(19) 

x,y= 0  x=0  y=0 

We  simplify  the  notation  by  defining  {/,m;=0;x  =  0,-,m-1}  ={M_,  =I,{Mi>0}  =  {M_l,M0,Ml,...}} 

=  {Ma,{_l  ia0j}  (a  set  of  m+ 1  linear  independent  matrices)  with  duals 
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satisfying  the  trace  orthogonality  conditions  Tr  [MaM p]=8a  p,  and  similarly  for 

{l,{Mf0}}  ->  {Np= Hj,0!}.We  therefore  write  (19)  as 


m— 1  m— 1  m— 1 

0=Yj  ®/  +  X  P7*’°^  (20) 

!,/=0  i=0  7=0 

using  the  abbreviations  pff  =  P(a  =  0, b  =  0 1 x  =  i, y  =  j),  if-’  =  p(a  =  0\x  =  i) and  P‘-°  =  p(b  =  0 1 y  =  j). 
For  the  measurement  matrices  M_x=i  =  i2x2,  and  Mii0  =\/2{i+mi  •  <?),  |w,|  <  1,  the  dual  matrices  are 
given  explicitly  by  M_x  =1  =  1/2(7-^  m,.  -d)  =  \/2(I -m-d),  and  M;  =  m.  ■  d  where  ni  ■  nr  =  Sij ,  \ni |  >  1, 

with  the  orthogonality  relations  Tr[l]  =  1,  Tr[M j  ]  =  o,  Tr[Mi  /]  =  0,  and  Tr[MiMj  ]  =  Sr  Using  the 
relationship  Tr[X  ®  Y]  =  Tr[X]  Tr[Y]  it  is  straightforward  to  verify  that  Tr[0]=l  and,  for  example, 
p"f  =Tr[OMi  ®  Xi  ]  which  picks  out  the  term  Mt  ®  /V  in  (20).  Other  probabilities  are  obtained  for 
example  as  Pfj  =  Tr[OMi  0(1- Nj )]  =  Tr[OM ,  01]- Tr[OMj  0Nj]=  Pf’  - Pff  =  £  iff  ~ iff 
=  p‘f  =  p(a  =  o,b  =  \\x  =  i,y  =  j).  Substituting  the  explicit  expressions  for  the  dual  matrices  into  (20) 
yields  the  general  expression  for  O  in  terms  of  products  of  Pauli  matrices 

1  n.  \ 

0  =  -[  2<®'  +  Pp)  + 1) (m,  ■  a ) ® % .  •  <x) 

,J=0  (21) 

m- 1  m- 1  v  J 

+  Z(2  X'  - 1)  (mt  •  <r)  0  7  +  -1)  I ®{nj  •  cr)  +  7  0  7  ]. 

i= 0  7=0 

Specializing  to  the  PR  correlations  in  (10)  given  by  P( «, />  |  x  =  /, y  =  j)  =1/2  <5^  i  /m()d2 
^Pff=\/2S0J.jmoil  with  marginals  pf’  =  Pf0  =1/2  v/,/,  yields  the  expression  for  the  NSPR  operator 


°pr  =^[  (me-a;)0(ne-<T)  +  (me-<T)0(no-(T)  +  (OTo-<T)0(ne-a;)-(mo-a;)0(no-<T)  +  7 07  ], 
where  me  =  £  m2i,  mo  =  £  m2/+1,  ne=  Y  &  =  Z  52y+i- 


(22) 


(=0,1,2,... 


(=0,1,2,... 


2=0, 1,2,... 


7=0, 1,2,... 


In  (22)  the  subscripts  {e,o}  denote  {even,odd}for  the  summation  over  even  and  odd  dual 
measurement  vectors.  Note  that  in  (22)  the  “single-c”  terms  <7.®/ and  / ® <7 .  (representing 

measurements  by  Alice  or  Bob  alone,  respectively)  have  dropped  out  since  the  marginal 
distributions  P(a\x)=P(b\y)=l/2  are  independent  of  a,b,x,y.  This  leaves  only  the  solely  two-party 
correlation  terms  cr.  ®cr.  and  the  maximally  mixed  term(/®/)/4.  For  the  bipartite  case  n— 2  often 

considered  in  the  literature  for  two  qubits,  each  with  two  measurement  directions  xe{m0,m,}  for 
Alice  and  y  =  {«„,«,}  for  Bob  (i.e.  a,b,x,y  =  {0,1})  we  obtain  the  simplified  form 
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0'PR  =  ( m0  ■  <t)  ®  («0  •  <t)  +  (m0  ■  <t)  ®  ( nx  ■  d)  +  ( mx  •  cf)  ®  («0  •  a) - (ihx  ■  5) ®  (nx  ■  5)  + 1  ®  /  ].  (23) 

Using  the  procedure  for  calculating  probabilities  discussed  after  equation  (4)  ,  the  following 
probabilities  can  be  computed  from 


«o 

«o 

n 

:o  «l 

«o 

«0 

m0 

"1/2 

1/2] 

m0 

"1/2 

1/2 

m0 

“0  0  “ 

m0 

“0 

0  " 

m0 

“1  1  “ 

Wj 

.1/2 

0 

mx 

.1/2 

0 

777, 

_°  1/2. 

Wj 

0 

1/2  J 

mx 

1  -1 

P(a=0,  b=0  |  x=ini ,  y=n. ),  P(a  =  l,b  =  1 1  mnn.\  P(a  =  0,b  =  l\mi,nj),  P(a  =  l,b  =  0  |  mi,nj),  E(mi,nj ) 


(24) 


a  ©  b  =  0, 

F  =  l/2for(x,j)G{(0,0),(0,l),(l,0)}, 


a  ©  b  - 1, 

P  =  l/2for(x,y)e{(l,l)}. 


Here,  the  correlations  in  (8)  are  computed  as  (see  (5)) 

E{m.,n.')i  ^.s{0 1}  =  P(a=Q,b=Q  \  x=m.,y=n.)  +  P(a  =  l,b  =  1 1  rh.,n.)  -  P(a  =  0  ,b  =  1 1  m.,n.)-  P(a  =  l,b  =  0 1  (25) 

with  corresponding  S  parameter  (see  (5)) 

S  =  E(m0,n0)  +  E(m0,nl)  +  E(ml,n0)-E(ml,nl)  =  4  =  SAM,  (26) 


achieving  the  algebraic  maximum  value  S m  =  4. 

For  the  case  of  two  qubits  with  m= 3  measurement  vectors  x  e  {m0 ,  mt , m2 }  for  Alice  and  y  =  {n0,nl,n2} 
for  Bob  (i.e.  a,b  =  {0,1},  with  x,y  =  {0,1,2})  we  obtain  from  (22)  the  probabilities  and 
correlations 


a, 

/72 

n0 

«i 

n2 

n 

!0 

«2 

«o 

A 

«2 

«o 

A 

«2 

m0 

"1/2 

1/2 

1/2" 

m0 

"1/2 

1/2 

1/2" 

m0 

“0 

0 

0" 

m0 

“0 

0 

0" 

m0 

“1 

1 

f 

1/2 

0 

1/2 

,  mx 

1/2 

0 

1/2 

,  "h 

0 

1/2 

0 

,  mt 

0 

1/2 

0 

,  m. 

1 

-1 

1 

m2 

.1/2 

1/2 

1/2. 

m2 

.1/2 

1/2 

1/2. 

m2 

0 

0 

0 

m2 

0 

0 

0 

_ 

m2 

1 

1 

1 

P(a=0, b=0  |  x=mi , y=n. ),  P(a  =  1, b  -  1 1  m. fij ),  P(a  =  0,Z?  =  0  |  m.,n.),  P(a  -  l,b  -  0  |  mi,nj),  E(mi,nj ) 
a®b  =  0,  a  ©  b  =  1, 

=^>  P  =  1/2 for (x, y)  e  {(e, e), (e, o), (o, e) } ,  ^P  =  l/2 for (x,  v)  g  |(o,o)}. 

In  (27)  e  =  {0,2}  denotes  even  indices  of  the  measurement  directions  while  o  =  { 1 }  denotes  odd 
indices.  We  achieve  the  algebraic  maximum  for  the  S  parameter,  generalizing  (26)  defined  as 

S  =  E(me,ne)  +  E(me,n0)  +  E(m0,ne)~  E(m0,n0)  =  4  =  SAM.  (28) 
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Note  that  the  dimension  of  the  measurement  vectors  m.  is  set  by  the  dimension D  =  d2-\  of  the 
Hilbert  space  of  the  observer,  which  simply  states  that  any  (Z)  +  l)x(Z)+l)  matrix  can  be  written  in 
term  of  the  (£>  +  l)x(£>  +  l)  identity  matrix  and  the  D  generators  of  su(d).  For  qubits,  D= 3  and  the 
three  generators  of  su( 2)  are  the  usual  Pauli  matrices  a.  For  a  given  set  of  in  measurement  3- 
vectors  {w,}  (vectors  in  the  Bloch  sphere,  \mt  |<l)  one  needs  to  solve  for  the  correspond  dual 
measurement  vectors  {mj }  satisfying  m,  • mj  =Sij.  We  write  these  equations  as  the  matrix  equation 
Mmx3  M3xm  =  imxm  where  the  ith  row  (i  =  {0,1,.. 1 })  of  (the  known  coefficient  matrix)  Mmx3 is 
and  the  yth  column  of  (unknowns)  M3xm  is  mr  By  linear  algebra,  there  exists  a  right  inverse  of 
Mmx3via  MRlghUm  =Mr(MMr)"‘  (if  (MMr)_1  exists)  if  the  columns  of  Mmx3  span  Rr" ,  which  can  only 

occur  for  m<D= 3.  The  systems  of  equations  is  under-determined  and  there  exists  at  least  one 
solution  (typically  and  infinite  number  due  to  undetermined  free  parameters).  This  is  the 
situation  for  probabilities  and  correlations  shown  in  (24)  and  (27)  for  the  case  m= 2  and  m= 3 
measurement  vectors,  respectively.  For  the  m>D= 3,  there  exists  at  most  one,  unique  solution  (if 
any).  This  is  the  least  squares  (LS)  solution  using  the  pseudo-inverse  Mmx3  given  by 
M/s  =  (MrM)"1Mr (if  (M'Mr  exits).  In  general,  the  LS  solution  has  non-zero  residual  errors  given 
by  Err=Mmx3MiS(3xm)  -lmxm,  corresponding  to  joint  probabilities  that  may  be  negative  for  some 

measurements  but  still  satisfy  the  (total  probability)  normalization  condition 
P(a,b\x,y)  =  \,  Vx, v.  Nonetheless,  it  is  instructive  to  perform  numerical  searches  in  the  case 

of  wj>3  of  random  measurement  vectors  to  seek  solutions  which  yield  all  joint  probabilities  in  the 
range  0<P(a,b\x,  v)<U  for  all  pairs  of  measurement  vectors  mnni  for  Alice  and  Bob  that  still  yield 

supra-correlations,  i.e.  0  <  S-Sq <  4-2^2  =  1.172. 

For  the  case  m= 4,  a  particular  solution  is  shown  in  (29)  that  yields  S-Sq=  0.102  (for  brevity,  we 
only  show  P(a=o, b=o  \  x=mi , y>=ni )  and  the  correlations In  general,  the  even/odd  structure  of 

the  correlations  E(m.,n.) 


n0 

«i 

n2 

«3 

n0 

«i 

n2 

«3 

m0 

'0.237 

0.395 

0.072 

0.406' 

m0 

'-0.052 

0.581 

-0.710 

0.623  ' 

mx 

0.162 

0.018 

0.381 

0.004 

fhx 

-0.350 

-0.927 

0.522 

-0.982 

m2 

0.469 

0.449 

0.341 

0.457 

m2 

0.875 

0.796 

0.365 

0.829 

m3 

0.249 

0.038 

0.481 

0.024 

m3 

-0.004 

-0.847 

0.923 

-0.905_ 

P(a=Q,  b= 0  |  x=ih! ,  y=Hj ),  E{mj,nj ) 


exhibited  in  the  cases  m< 3  ((24)  and  (27))  is  destroyed,  yet  they  still  produce  supra-correlations 
S-Sq  >0.  For  each  value  of  m  in  Figure  24  (left)  we  searched  105  random  trials  of  the 


Approved  for  Public  Release;  Distribution  Unlimited. 

34 


measurement  vectors  {m,.,ny.}Ue{0ili  >m_1}  and  plot  the  value  of  S-Sq  for  the  first  solution  encountered 
in  which  (i)  we  find  proper  joint  probability  distributions  0<P(a, b  \  x=mi , y=nj  )<  1  for  all 

measurement  vectors,  and  (ii)  which  produce  supra-correlations,  S-Sq> 0.  In  Figure  24  (middle), 
we  plot  the  minimum 


S-SQ 


Figure  24.  Numerical  simulations  for  m={2,3,4,. .  .,12} measurement  vectors. 

eigenvalue  Amm  of  the  matrix  O  in  (22).  The  negative  value  of  Amin  indicates  that  O  is  not  realized 
by  a  proper  quantum  state  (i.e.  a  positive,  Hermitian  operator,  p>0).  The  rightmost  plot  in  Figure 
24  is  the  iteration  number  at  which  the  first  set  of  measurement  vectors  was  found  which 
produced  supra-correlations.  For  the  values  of  13  <  m  <  20  numerically  explored,  no  supra- 
correlations  solutions  were  found  within  105  trials  (the  plot  indicates  that  it  becomes 
exponentially  hard  to  find  such  a  solution). 


No  Signaling  (NS)  Correlations:  3-Qubits 

The  bipartite  results  of  the  previous  section  for  «=2-qubits  are  straightforwardly  extended  to  the 
tripartite  case  of  «=3-qubits  with  similar  implications.  Here  the  generalization  of  the  bipartite 
CHSH  nonlocality  parameter  S  is  given  by  the  Svetlichny  [Svetlichny87]  inequality  (SI)  relating 
correlations  E(a,b,c\x,y,c)  between  three  spacelike  separated  parties  A,  B,  C 

S  =E(a,b,c  |  0,0,0 )  +  E(a,b,c  |  0,1, 0)  +  E(a,b,c\  1,0, 0)  -E{a,b,c  1 1,1,0) 

+  E(a,b,c  1 0,0,1  )-E(a,b,c  |  0,1,1  )-E{a,b,c  1 1,0,1  )-E(a,b,c  \  1,1,1).  ^  ^ 

The  SI  has  the  bounds  (i)  \s\  <  Sc  =  4  for  classical  correlations,  (ii)  „s)  <SQ=  4V2  for  quantum 
correlations,  with  (iii)  the  algebraic  upper  bound  given  by  |.sj  <  SAM  =  8,  achieved  when  the 

correlations  in  (14)  take  the  values  E=  1  if  they  are  preceded  by  a  plus  sign,  and  E=- 1  if  they  are 
preceded  by  a  minus  sign.  The  generalization  of  the  PR  correlations  of  (30)  is  given  by 
[Xiangll] 


TPR  Box: 


P(a,b,c  |  x,y,  z )  = 


if  a®b®c  =  x-y®y-z®x-z 
otherwise 


(31) 
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often  referred  to  as  a  tripartite  PR  (TPR)  box.  The  marginal  distributions  of  (15)  are  again 
isotropic  and  satisfy  the  NS  constraint,  i.e.  P(a,b\x,y,z )  =  P(a, b\x,y )=l/4  for  all  a,b,x,y,z  and 
P(a,  \x,y)=P(a\x)=l/2  for  all  a,x,y,  and  similarly  for  all  other  marginal  probability  distributions. 

For  the  case  of  n= 3  qubits  (r= 2  output  measurement  values)  a,b,c  ={0,1),  with  m  possible 
measurement  vectors  for  each  observer,  x,y,z  ={0,l,...,m-l)  we  again  find  that  only  the  highest 
(three  party)  correlations  term  and  the  maximally  mixed  term  are  non-zero  in  the  expression  for 
Otpr 

Otpr  =  ]-[  {(jhe-G)®(ne-a)  +  (me-a)®(n0-a)  +  (jh0-a)®(ne-G)-(jh0-a)®(n0-a)}®{re-a) 

o 

~{{m0  ■a)®(n0-a)  +  (m0-a)®(ne-a)  +  (me-a)®(n0  ■a)-(me-a)®(ne-a)}®(ro  •<r)  +  /®/®/],  (32) 
where  qe  =  £  q2i,  qo=  £  q2M,  q  =  {m,n,r}. 

i= 0,1,2,...  ~  i=0,l,2,... 

The  regular,  even/odd  (mod  2)  structure  of  Otpr  in  (32)  reflects  the  non-zero  structure  of  the  TRP 
probabilities  in  (31),  and  can  be  seen  as  an  additional  single  qubit  generalization  of  0PR  in  (22). 
That  is,  the  2-qubit  term  in  the  first  curly  brackets  in  (32) 
{( rne  ■  a)  ®  ( ne  •  <?)  +  ( rne  ■  cr)  ®  ( n0  •  <?)  +  ( rh0  •  <?)  ®  ( ne  •  <?)  -  (mB  ■  cr)  ®  ( n0  ■  <?)}  tensor-producted  with  the 
remaining  “even”  qubit  term  (?e  •  a) ,  is  precisely  two-party  correlation  term  that  appears  in  0PR 
in  (22).  Similarly,  the  term  in  the  second  curly  bracket  in  (32)  {(m0  ■  a)  ®  (n0  ■  a)  +  (m0  ■  a)  ®  (ne  ■  a) 
+(me-a)®(n0-a)-(me-a)®(ne-a)}  tensor-producted  with  the  remaining  “odd”  qubit  term  (Fo-<?) 

(with  the  accompanying  minus  sign)  is  just  the  bit  flip  ( e<->o )  of  the  previous  two-party 
correlation  term.  Again,  we  can  achieve  the  algebraic  maximum  Sam=  8  when  each  party  has  (for 
the  case  of  qubits)  at  most  m= 3  measurement  vectors  (for  exactly  the  same  linear  algebraic 
reason  for  the  n= 2  bipartite  case).  Further,  as  in  the  bipartite  case,  we  can  find  particular  NS 
supra-correlation  solutions  0<.S-.S);<4-2V2  for  m> 3,  but  which  become  increasingly  hard  to  find 
the  larger  the  value  of  m. 

4.6  Ongoing  cluster  state  algorithm  research 

As  part  of  our  theoretical  research  goals  for  this  project,  we  have  been  investigating  quantum 
algorithm  development  in  the  cluster  state  paradigm.  In  almost  every  implementation  of  quantum 
computing  one  of  the  first  algorithms  that  people  have  developed  and  tested  is  an  unstructured 
quantum  search  algorithm.  In  the  photonic  circuit  model  of  quantum  computing  [DiVincenzoOO] 
this  algorithm  is  structurally  very  similar  to  that  first  proposed  by  Grover  and  as  such  is  called 
Grover's  algorithm  [Grover97].  Since  the  development  of  MBQC,  and  particularly  in  the 
photonic  implementation  of  MBQC,  also  called  the  cluster  state  model,  parallels  have  been 
drawn  between  the  photonic  circuit  model  and  MBQC  model.  This  includes  comparing  the 
circuit  model  of  Grover's  algorithm  to  a  cluster  state  model  of  a  four  element  quantum  search 
[Walther05],  Given  the  difficulty  of  creating  large  cluster  states  this  four  element  search  is  the 
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largest  that  has  been  performed  to  date.  The,  2"=22=4  element  search  requires  «2=22=4  qubits 
arranged  in  a  square  or  "box"  cluster  state  [Walther05],  The  next  larger  2'=8  element  search 
requires  a  significantly  more  difficult  to  produce  3x3  square  cluster  state  and  so  on. 

Using  the  4-qubit  box  state,  as  illustrated  in  Figure  25,  the  measurement  based  quantum  search 
appears  to  be  very  similar  to  that  of  Grover.  The  input  state  is  similar,  the  desired  output  state  is 


Figure  25.  Square  cluster  states  (a)  (n= 2)  4-element  “box,”  (b)  (n= 3)  9-element. 

similar,  the  number  of  iterations  is  the  same  and  the  oracle's  tagging  operation  appears  similar. 
This  has  led  many  people  to  call  the  MBQC  quantum  search,  Grover's  algorithm  [Grover97], 
However,  there  is  often  a  small  caveat  that  is  overlooked.  As  Zeilinger  et  al.  wrote  “Remarkably, 
the  inversion-about-the-mean  process  is  'hard-wired'  into  the  structure  of  the  cluster  state  and  is 
automatically  implemented”  [Walther05].  We  found  this  automatic  implementation  of  an 
inversion  about  a  mean,  or  equivalently  amplitude  amplification  (when  combined  with  tagging), 
to  be  an  imprecise  claim  and  as  such  have  investigated  the  measurement  based  quantum  search 
algorithm  numerically  for  larger  systems.  We  have  found  that  a  significant  number  of  non-trivial 
differences  exist  between  the  standard  description  of  Grover's  algorithm  and  its  implementation 
in  the  circuit  model  and  the  MBQC  search  algorithm  for  larger  systems.  This  leads  us  to 
question  whether  the  MBQC  search  algorithm  is  some  variant  of  Grover's  algorithm,  or  rather  a 
different  search  algorithm  all  together. 

Two  cluster-states  able  to  perform  an  unsorted  search  for  a)  the  (n= 2)  22=4  element  search  on 
Y^box)  ar|d  b)  the  (n= 3)  23=8  element  search  on  |C9)  are  shown  in  Figure  25.  Each  circle  is  a 
single  qubit  initialized  in  the  |+)  state.  The  solid  lines  connecting  the  circles  indicate  the  action 
of  a  CZ  gate  between  these  qubits,  i.e.  entanglement  between  the  qubits.  Logical  qubits  are 
arranged  in  rows  of  the  grid,  initially:  (1,3)  and  (2,4)  for  |  Vbox)’  anc*  (1,4,7),  (2,5,8),  and  (3,6,9) 
for  |C9).  Figure  25(a)  shows  the  trivial  MBQC  measurements  using  one  application  of  an  iterator 
Ii  (large  red  dashes  acting  on  columns  of  qubits)  on  the  4-qubit  box  state.  Figure  25(b)  Shows  the 
two  iterations  for  the  eight  element  search  I\,h.  All  of  the  measurements  in  an  iterator  /j  can  be 
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implemented  simultaneously.  The  outcomes  of  each  7j  determine  the  basis  of  the  next  iteration 
w  (Equivalently  7j  determines  the  Pauli  correction  factors  that  can  be  applied  latter  in  post¬ 
processing).  The  output  state  (small  green  dashes)  holds  the  logical  state  at  the  end  of  the 
algorithm.  The  effective  logical  qubits  after  the  first  iteration  I\  are  (3)  and  (4)  (the  output 
states)  for  | and  (4,7),  (5,8),  and  (6,9)  for  |C9). 


Measurement  Sequence  ( X\  ;  I2) 

Search  Element  i 

Output  State  I’J'out.t) 

(0,0,0  ;  0,0,0) 

1 

(1,0, 0,0, 0,0, 0,0) 

(0,0, 7T  ;  0,7T,0) 

2 

(0,1, 0,0, 0,0, 0,0) 

(0,0,0  ;  0,0, 7r) 

3 

(0,0,1,0,0,0,0,0) 

(0,0,7 r  ;  0,7T,7r) 

4 

(0,0,04,0,0,0,0) 

(0,0, 7T  ;  0,0,0) 

5 

(0,0,0,04,0,0,0) 

(0,0,0  ;  0,tt,0) 

6 

(0,0,0,0,04,0,0) 

(0,0, 7r  ;  0,0, 7r) 

7 

(0,0,0,0,0,04,0) 

(0,0,0  ;  0,7T,7T ) 

8 

(0,0,0,0,0,0,04) 

Figure  26.  Measurement  patterns  and  resulting  output  for  8-element  MBQC  search. 

Here  we  pause  to  note  something  rather  undesirable.  If  we  consider  the  measurement  of  the 
vertical  columns  to  be  a  “step”  in  the  iteration  (i.e.  1,2,3  or  4,5,6  from  Figure  25(b)),  then  the 
MBQC  algorithm  does  not  have  a  fixed,  constant  iterator.  The  iterator  will  in  general  vary 
between  steps  and  it  is  neither  sequential  nor  random  as  can  be  seen  from  the  first  column  of 
Figure  26.  The  various  iterators  must  be  carefully  chosen  by  the  oracle  at  each  step.  This  appears 
very  different  from  Grover's  iterator.  In  Grover's  algorithm  the  iterator  is  chosen  by  the  oracle  at 
the  start  of  the  calculation  and  then  remains  a  constant  regardless  of  the  number  of  iteration 
steps.  Figure  26  shows  measurement  settings  for  the  6-qubit  measurements  that  result  in  each  of 
the  eight  possible  outputs  of  the  MBQC  search.  Note  the  output  state  is  corrected  for  overall 
phase,  so  the  tags  are  not  unique.  In  addition,  all  measurements  presented  are  presumed  to  have 
given  the  “correct”  output  (i.e.  only  entries  in  section  1  of  YV0Ut)  are  populated),  thus  negating 
the  need  for  trivial  feed  forward  corrections  typical  in  MBQC.  We  have  found  a  way  of 
modifying  the  MBQC  search  algorithm  such  that  the  iterator  is  a  constant  not  only  between 
steps,  but  is  also  constant  for  any  desired  output.  This  work  an  ongoing  project  for  FY13,  and  the 
research  results  are  being  written  up  for  submission  to  Physical  Review  A. 

5.0  CONCLUSIONS 

Multipli-entangled  photons  from  a  spontaneous  parametric  down-conversion  source 

This  report  describes  research  on  the  Schioedtei  source,  a  unique  type  II  SPDC  source  design  for 
which  additional  in-depth  information  can  be  obtained  through  our  previously  published  papers 
[Fantoll,  Peters  12],  Schioedtei  generates  up  to  six  pairs  of  entangled  photons  per  pass  through 
the  type  II  crystal  assembly.  This  configuration  surpasses  the  typical  single  entangled  pair 
generated  per  pass  found  in  standard  type  II  SPDC  sources.  Concurrently  Schioedtei  generates  a 
variety  of  states  atypical  of  being  produced  from  a  single  photon  source.  Useable  photon 
generation  rates  (two  and  four  photon)  have  been  observed,  thus  showing  its  feasibility  as  a 
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direct  generation  source  of  entangled  photons  for  quantum  optics/entanglement  experiments. 
The  six  pairs  of  photons  produced  are  directly  applicable  to  the  generation  of  linear,  box, 
butterfly  and  a  multitude  of  other  cluster  states.  The  utility  of  the  Schioedtei  source  is  (i)  its 
reduced  experimental  footprint  compared  to  standard  multi-crystal/multi-pass  experiments,  (ii)  it 
generates  a  variety  of  entangled/separable  states,  and  (iii)  generated  states  are  amenable  towards 
cluster  state  generation.  Furthermore,  the  generated  photons  from  Schioedtei  are  the  input  states 
for  our  bulk  optical  gates  and  QPICs. 

A  multi-layer  three  dimensional  superconducting  nanowire  photon  detector 

The  multilayer  superconducting  number-resolving  photon  detector  represents  a  significant 
improvement  on  current  single  layer  meander  devices.  The  device  will  have  significantly  higher 
number  resolution,  while  maintaining  a  useful  detection  area.  It  has  several  parameters  which 
can  control  the  reset  time  to  avoid  latching  while  still  minimizing  the  rest  time.  An  array  of 
pixels  of  arbitrary  number,  size  and  shape  is  possible.  The  active  area  of  the  detector  can  be 
tuned  by  changing  the  number  or  the  shapes  of  the  pixels.  The  fill  factor  of  the  detector  should 
be  at  least  equal  to  that  of  current  nanowire  meanders  and  given  the  potential  reduction  of  the 
current  crowding  effect  significantly  higher.  As  a  final  note  we  will  point  out  that  the  multi-layer 
superconducting  number-resolving  photon  detector  can  also  give  a  rough  spatial  distribution  of 
the  incident  photons.  These  advantages  are  compelling  evidence  for  the  construction  and  testing 
of  multi-layer  superconducting  number-resolving  photon  detectors. 

Theory/experimental  requirements  of  imperfect  two-qubit  linear  optical  photonic  gates 

We  have  shown  the  theoretical  basis  and  interest  for  this  experiment.  At  this  time  it  is  the  only 
apparent  means  of  experimentally  confirming  the  numerical  data  presented  above,  which 
quantifies  the  trade-off  between  fidelity  and  success,  for  the  CZ  or  CNOT  gate.  The  experimental 
setup  may  naturally  be  extended  to  explore  the  behavior  of  other  quantum  gates  of  interest.  The 
components  needed  for  the  execution  of  the  experiment  are  well  within  the  means  of  many 
experimental  groups.  The  main  stumbling  block  is  the  expense  of  purchasing  number-resolving 
detectors.  However,  any  group  already  possessing  these  detectors  should  be  able  to  implement 
this  scheme  with  relative  ease. 

Nonlocality,  entanglement  witnesses  and  supra-correlations 

In  this  area  of  research  we  have  examined  the  structure  of  supra-correlations  that  are  stronger 
than  quantum  and  hence  not  realizable  by  a  physical  (positive)  quantum  state  p>0.  The  supra- 
correlations  are  intriguing  because  they  arise  from  valid  probability  distributions,  first  put  forth 
by  Popescu  and  Rohrlich  (PR),  that  satisfy  the  no-signaling  principle  of  special  relativity  as  well 
as  all  the  usual  normalization  condition  on  the  joint  and  marginal  distributions.  Thus,  the  fact  that 
nature  is  not  able  to  realize  these  supra-correlations  points  to  hidden  structure  underlying  how 
quantum  correlations  can  be  distributed  amongst  spacelike  separated  parties.  Our  work  has 
examined  the  structure  and  distribution  of  PR  correlations  in  2-  and  3 -qubit  systems  by  explicitly 
constructing  “states”  (not  necessarily  positive  quantum  states)  that  exhibit  supra-correlations  for 
a  fixed,  but  arbitrary  number,  of  measurements  available  to  each  party.  We  have  shown  that  the 
PR  correlations  involve  only  solely  n-party  correlations  amongst  the  n  observers.  We  have 
extended  this  study  to  include  n-party  correlations  that  capture  the  essential  features  of  the  PR 
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correlations  and  do  not  rely  on  predetermined  measurements  between  the  n  participants. 
Additionally,  by  constructing  constraints  based  on  the  positivity  and  purity  of  an  arbitrary  n-qubit 
state  we  have  shown  the  “unreasonableness”  of  the  PR  correlations  in  that  they  encode  more 
correlations  than  are  physically  allowed  by  nature  [see  details  in  Alsingl2],  In  future  work  we 
will  couple  this  approach  of  studying  how  correlations  are  distributed  amongst  the  n  parties  to 
the  study  of  quantum  entanglement.  The  study  of  entanglement  [Horodecki09]  is  an  important, 
but  difficult  field,  only  well  understood  for  the  case  of  two  qubits  (both  pure  and  mixed),  and  to  a 
lesser  degree,  for  pure  3-qubit  systems.  A  fruitful  area  to  investigate  next  are  pure  3-qubit 
systems,  where  a  generalized  (though  non-unique)  Schmidt  decomposition  holds  [AcinOO].  We 
purport  that  an  examination  of  the  distribution  of  correlations,  bounded  by  physically  imposed 
constraints  on  e.g.  positivity  and  purity,  coupled  with  the  description  of  entanglement  in  terms  of 
the  tangle,  as  initiated  in  this  work,  can  shed  further  light  on  the  classification  of  pure  tripartite 
systems. 
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7.0  LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


2D 

2  Dimensional 

3D 

3  Dimensional 

a-BBO 

Alpha  barium  borate 

AFRL 

Air  Force  Research  laboratory 

APD 

Avalanche  photodiode 

a-WxSii.x 

Amorphous  Tungsten  Silicon 

B-BBO 

Beta  barium  borate 

BiBO 

Bismuth  borate 

BI 

Bell  Inequality 

CCD 

Charged  coupled  device 

CCM 

Coincidence  counting  module 

CHSH 

Clauser-Home-Shimony-Holt 

CJI 

Choi-Jamiolkowski  isomorphism 

CNOT 

Controlled  NOT  (gate) 

CW 

Continuous  wave 

CZ 

Controlled  Z  (gate) 

fs 

femtosecond 

GVM 

Group  velocity  matching 

HP 

High  power 

JAG 

Judge  Advocate  General 

LS 

Least  squares 

MBQC 

Measurement  based  quantum  computing 

MgO 

Magnesium  Oxide 

MHz 

megahertz 

mm 

millimeter 
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NbN 

Niobium  Nitride 

NbTiN 

Niobium  Titanium  Nitride 

NIST 

National  Institute  for  Standards  and  Technology 

nm 

nanometer 

NS 

No-signaling 

0P0 

Optical  Parametric  Oscillator 

PQNS 

Post  Quantum  No-signaling 

PR 

Popescu-Rohrlich 

QM 

Quantum  Mechanics 

QIS 

Quantum  Information  Science 

QPIC 

Quantum  Photonic  Integrated  Circuit 

SHG 

Second  harmonic  generation 

Si 

Silicon 

Si-APD 

Silicon  avalanche  photodiode 

SI 

Svetlichny  Inequality 

SNSPD 

Superconducting  Nanowire  Single  Photon  Detector 

SPCM 

Single  photon  counting  module 

SPDC 

Spontaneous  parametric  downconversion 

SR 

Special  Relativity 

THG 

Third  harmonic  generation 

TPR 

Tripartite  Popescu-Rohrlich 

UV 

Ultraviolet 

Approved  for  Public  Release;  Distribution  Unlimited. 

46 


